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International Filing Date 
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Title of Invention 

PARTICLE ANALYSIS SYSTEM AND 
METHOD 



Authorized Officer: 
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Mailing Date 
28 January 1998 

Agent's File 

Reference: 
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RESPONSE TO INVITATION TO CORRECT DEFECTS 



Commissioner of Patents and Trademarks 
Box PCT 

Washington, D.C. 20231 

Attn: RO/US 

Dear Sir/Madam: 

In response to the Invitation to Correct Defects in the International 
application mailed 29 December 1997, Applicant submits herewith eight 
duly executed Appointments of Agent, new substitute description, claims 
and abstract and Sheet 3 of the Request correcting the spelling of Jeffrey 
Kao's first name. In the course of correcting the Description, a printing 
error was noted regarding equations (10) and (18). Specifically, stray 
characters were printed instead of the intended ellipses. These errors 



"Express WtafT tofcel £> frvt ~> C^lk ^ 5 A.^Q/l/^ . 

Date of Oepostt ^ ^ /U^^^^ / X 

I hereby certify thai this paper or fee/ is being deposited with 
the United States Posts* Service "Express Mai! Pest Office to 
' ^dressoe' service under 37CFR § 1.10 on the date indicated 
• 9 end is addressed to 0)6 Assistant Corj^missioner for Patents, 
Arlington, D.C. 2C231 J 




As 



Signature or parson n*ail;#Lj paper or fee 



r 



C 




have also been corrected. The ellipses were properly included in the same 
equations of the U.S. Patent application Serial No. 08/747,112 from which 
priority is claimed and which was incorporated by reference in the PCT 
application as originally filed on page 56 and as amended on page 64. 
Thus it is believed that this subject matter is not being newly added. We 
believe that this application conforms to the requirements of the PCT. 



Respectfully submitted 




L. Scott Paynter, #39,797 
Woodard, Emhardt, Naughton, 



Moriarty & McNett 
Bank One Center/Tower, Suite 3700 
111 Monument Circle 
Indianapolis, Indiana 46204 US 
(317) 634-3456 



Enclosures: Eight Appts. Of Agent, substitute 

pages of Description, Claims and Abstract and 
Sheet 3 of Request 
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PATENT COOPERATION TREATY 



Appointment of Agent or Common Representative 



The undersigned applicant, PURDUE RESEARCH FOUNDATION, hereby appoints 
Harold R. Woodard, No. 16,214; C. David Emhardt, No. 18,483; Joseph A. Naughton, Jr., 
No. 19,814; John V. Moriarty, No. 26,207; John C. McNett, No. 25,533; Thomas Q. Henry, 
No. 28,309; James M. Durlacher, No. 28,840; Charles R. Reeves, No. 28,750; Vincent O. 
Wagner, No. 29,596; Steve Zlatos, No. 30,123; Spiro Bereveskos, No. 30,821; William F. 
Bahret, No. 31,087; Clifford W. Browning, No. 32,201; R. Randall Frisk, No. 32,221; 
Daniel J. Lueders, No. 32,581; Michael D. Beck, No. 32,722, Kenneth A. Gandy, No. 
33,386; Timothy N. Thomas, No. 35,714; Kerry P. Sisselman, No. 37,237; Kurt N. Jones, 
No. 37,996; John H. Allie, No. 39,088; Jeffrey A. Michael, No. 37,394; Holiday W. Banta, 
No. 40,311; Troy J. Cole, No. 35,102; L. Scott Paynter, No. 39,797; J. Andrew Lowes, No. 
40,706; Darrin Wesley Harris, No. 40,636; Matthew R. Schantz, No. 40,800; Gregory B. 
Coy, 40,967; Lisa A. Hiday, No. 40,036; John V. Daniluck, 40,581 and James L. Rowe, 
No. 18,448, all of the same address and who are members of the Bar of the State of 
Indiana, constituting the firm of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, Indianapolis, Indiana 46204 United 
States of America 



to act on its behalf before the competent International Authorities in connection with the 
international application concerning PARTICLE ANALYSIS SYSTEM AND METHOD, 
Agent's File Reference PUR-63/75P, International Application Number unknown at this 
time, filed with the United States receiving office and to receive payments on its behalf. 
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Secretary 
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PATENT COOPERATION TREATY 
Appointment of Agent or Common Representative 



The undersigned applicant, Eva SEVICK-MURACA, hereby appoints Harold R. 
Woodard, No. 16,214; C. David Emhardt, No. 18,483; Joseph A. Naughton, Jr., No. 
19,814; John V. Moriarty, No. 26,207; John C. McNett, No. 25,533; Thomas Q. Henry, No. 
28,309; James M. Durlacher, No. 28,840; Charles R. Reeves, No. 28,750; Vincent O. 
Wagner, No. 29,596; Steve Zlatos, No. 30,123; Spiro Bereveskos, No. 30,821; William F. 
Bahret, No. 31,087; Clifford W. Browning, No. 32,201; R. Randall Frisk, No. 32,221; 
Daniel J. Lueders, No. 32,581; Michael D. Beck, No. 32,722, Kenneth A. Gandy, No. 
33,386; Timothy N. Thomas. No. 35,714; Kerry P. Sisselman, No. 37,237; Kurt N. Jones, 
No. 37,996; John H. Allie, No. 39,088; Jeffrey A. Michael, No. 37,394; Holiday W. Banta. 
No. 40,311; Troy J. Cole, No. 35,102; L. Scott Paynter, No. 39,797; J. Andrew Lowes, No. 
40,706; Darrin Wesley Harris, No. 40,636; Matthew R. Schantz, No. 40,800; Gregory B. 
Coy, 40,967; Lisa A. Hiday, No. 40,036; John V. Daniluck, 40,581 and James L. Rowe, 
No. 18,448, all of the same address and who are members of the Bar of the State of 
Indiana, constituting the firm of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, Indianapolis, Indiana 46204 United 
States of America 



to act on his/her behalf before the competent International Authorities in connection with 
the international application concerning PARTICLE ANALYSIS SYSTEM AND METHOD, 
Agent's File Reference PUR-63/75P, International application Number unknown at this 
time, filed with the United States receiving office and to receive payments on his/her 
behalf. 

(Place) (Meri Lej^yttu^ J^a^ta^cJ (Date) /O^m^^o* 3 , 



[X] as agent [ ] as common representative 





Name: Eva SEVICK-MURACA 



PATENT COOPERATION TREATY 
Appointment of Agent or Common Representative 

The undersigned applicant, Joseph PIERCE, hereby appoints Harold R. Woodard, 
No. 16,214; C. David Emhardt, No. 18,483; Joseph A. Naughton, Jr., No. 19,814; John V. 
Moriarty, No. 26,207; John C. McNett, No. 25,533; Thomas Q. Henry, No. 28,309; James 
M. Durlacher, No. 28,840; Charles R. Reeves, No. 28,750; Vincent O. Wagner, No. 
29,596; Steve Zlatos, No. 30,123; Spiro Bereveskos, No. 30.821 ; William F. Bahret, No. 
31,087; Clifford W. Browning, No. 32,201; R. Randall Frisk, No. 32,221; Daniel J. 
Lueders, No. 32,581; Michael D. Beck, No. 32,722, Kenneth A. Gandy, No. 33,386; 
Timothy N. Thomas, No. 35,714; Kerry P. Sisselman, No. 37,237; Kurt N. Jones, No. 
37,996; John H. Allie, No. 39,088; Jeffrey A. Michael, No. 37.394; Holiday W. Banta, No. 
40,311; Troy J. Cole, No. 35,102; L. Scott Paynter, No. 39,797; J. Andrew Lowes, No. 
40,706; Darrin Wesley Harris, No. 40,636; Matthew R. Schantz, No. 40,800; Gregory B. 
Coy, 40.967; Lisa A. Hiday, No. 40.036; John V. Daniluck. 40,581 and James L. Rowe, 
No. 18,448, all of the same address and who are members of the Bar of the State of 
Indiana, constituting the firm of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, Indianapolis, Indiana 46204 United 
States of America 

[X] as agent [ ] as common representative 

to act on his/her behalf before the competent International Authorities in connection with 
the international application concerning PARTICLE ANALYSIS SYSTEM AND METHOD, 
Agent's File Reference PUR-63/75P, International application Number unknown at this 
time, filed with the United States receiving office and to receive payments on his/her 
behalf. 

(Place) J^^b^l^ uJ>J^ (Date) l^/l^/^l- 
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PATENT COOPERATION TREATY 
Appointment of Agent or Common Representative 



The undersigned applicant, Steven RICHTER, hereby appoints Harold R. Woodard. 
No. 16.214; C. David Emhardt, No. 18,483; Joseph A. Naughton, Jr., No. 19,814; John V. 
Moriarty, No. 26,207; John C. McNett, No. 25,533; Thomas Q. Henry, No. 28,309; James 
M. Durlacher, No. 28,840; Charles R. Reeves, No. 28,750; Vincent O. Wagner, No. 
29,596; Steve Zlatos, No. 30,123; Spiro Bereveskos, No. 30,821; William F. Bahret, No. 
31,087; Clifford W. Browning, No. 32,201; R. Randall Frisk, No. 32,221; Daniel J. 
Lueders, No. 32,581; Michael D. Beck, No. 32,722, Kenneth A. Gandy, No. 33,386; 
Timothy N. Thomas, No. 35,714; Kerry P. Sisselman, No. 37,237; Kurt N. Jones, No. 
37,996; John H. Allie, No. 39,088; Jeffrey A. Michael, No. 37,394; Holiday W. Banta, No. 
40,311; Troy J. Cole, No. 35,102; L. Scott Paynter, No. 39,797; J. Andrew Lowes, No. 
40,706; Darrin Wesley Harris, No. 40,636; Matthew R. Schantz, No. 40,800; Gregory B. 
Coy, 40,967; Lisa A. Hiday, No. 40,036; John V. Daniluck, 40,581 and James L. Rowe, 
No. 18,448, all of the same address and who are members of the Bar of the State of 
Indiana, constituting the firm of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, Indianapolis, Indiana 46204 United 
States of America 



to act on his/her behalf before the competent International Authorities in connection with 
the international application concerning PARTICLE ANALYSIS SYSTEM AND METHOD, 
Agent's File Reference PUR-63/75P, International application Number unknown at this 
time, filed with the United States receiving office and to receive payments on his/her 
behalf. 



[X] as agent [ ] as common representative 
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PATENT COOPERATION TREATY 
Appointment of Agent or Common Representative 

The undersigned applicant, Rajesh SHINDE, hereby appoints Harold R. Woodard, 
No. 16,214; C. David Emhardt, No. 18,483; Joseph A. Naughton, Jr., No. 19,814; John V. 
Moriarty, No. 26,207; John C. McNett, No. 25,533; Thomas Q. Henry, No. 28,309; James 
M. Durlacher, No. 28,840; Charles R. Reeves, No. 28,750; Vincent O. Wagner. No. 
29,596; Steve Zlatos, No. 30,123; Spiro Bereveskos, No. 30,821 ; William F. Bahret, No. 
31,087; Clifford W. Browning, No. 32.201; R. Randall Frisk, No. 32,221; Daniel J. 
Lueders, No. 32.581; Michael D. Beck, No. 32,722, Kenneth A. Gandy, No. 33,386; 
Timothy N. Thomas, No. 35,714; Kerry P. Sisselman, No. 37,237; Kurt N. Jones, No. 
37,996; John H. Allie, No. 39,088; Jeffrey A. Michael, No. 37,394; Holiday W. Banta, No. 
40,311; Troy J. Cole, No. 35,102; L. Scott Paynter, No. 39,797; J. Andrew Lowes, No. 
40,706; Darrin Wesley Harris, No. 40,636; Matthew R. Schantz, No. 40,800; Gregory B. 
Coy. 40,967; Lisa A. Hiday, No. 40,036; John V. Daniluck, 40,581 and James L. Rowe, 
No. 18,448, all of the same address and who are members of the Bar of the State of 
.Indiana, constituting the firm of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, Indianapolis, Indiana 46204 United 
States of America 

[X] as agent [ ] as common representative 

to act on his/her behalf before the competent International Authorities in connection with 
the international application concerning PARTICLE ANALYSIS SYSTEM AND METHOD, 
Agent's File Reference PUR-63/75P, International application Number unknown at this 
time, filed with the United States receiving office and to receive payments on his/her 
behalf. 

(Place) L'fc)p/A' bha^Jif IN (Date) ^^#5*^^ oil °l °\<? 
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PATENT COOPERATION TREATY 
Appointment of Agent or Common Representative 

The undersigned applicant, Ganesh BALGI, hereby appoints Harold R. Woodard, 
No. 16,214; C. David Emhardt, No. 18,483; Joseph A. Naughton, Jr., No. 19,814; John V. 
Moriarty, No. 26,207; John C. McNett, No. 25,533; Thomas Q. Henry, No. 28,309; James 
M. Durlacher, No. 28,840; Charles R. Reeves, No. 28,750; Vincent O. Wagner, No. 
29,596; Steve Zlatos, No. 30,123; Spiro Bereveskos, No. 30,821 ; William F. Bahret, No. 
31,087; Clifford W. Browning, No. 32,201; R. Randall Frisk, No. 32,221; Daniel J. 
Lueders, No. 32,581; Michael D. Beck, No. 32.722, Kenneth A. Gandy, No. 33,386; 
Timothy N. Thomas, No. 35,714; Kerry P. Sisselman, No. 37,237; Kurt N. Jones, No. 
37,996; John H. Allie, No. 39,088; Jeffrey A. Michael, No. 37,394; Holiday W. Banta, No. 
40,311; Troy J. Cole, No. 35,102; L. Scott Paynter, No. 39,797; J. Andrew Lowes, No. 
40,706; Darrin Wesley Harris, No. 40,636; Matthew R. Schantz, No. 40,800; Gregory B. 
Coy, 40,967; Lisa A. Hiday, No. 40,036; John V. Daniluck, 40,581 and James L. Rowe. 
No. 18,448, all of the same address and who are members of the Bar of the State of 
Indiana, constituting the firm of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, Indianapolis, Indiana 46204 United 
States of America 

[X] as agent [ ] as common representative 

to act on his/her behalf before the competent International Authorities in connection with 
the international application concerning PARTICLE ANALYSIS SYSTEM AND METHOD, 
Agent's File Reference PUR-63/75P, International application Number unknown at this 
time, filed with the United States receiving office and to receive payments on his/her 
behalf. A * ^ V=^7 

(Place) ^ W \ (Date) ^ 
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Name: Ganesh BALGI 




PATENT COOPERATION TREATY 
Appointment of Agent or Common Representative 

The undersigned applicant, Jeffery KAO, hereby appoints Harold R. Woodard, No. 
16,214; C. David Emhardt, No. 18,483; Joseph A. Naughton. Jr., No. 19,814; John V. 
Moriarty, No. 26,207; John C. McNett, No. 25,533; Thomas Q. Henry, No. 28,309; James 
M. Durlacher, No. 28,840; Charles R. Reeves, No. 28,750; Vincent O. Wagner, No. 
29,596; Steve Zlatos, No. 30,123; Spiro Bereveskos, No. 30,821; William F. Bahret, No. 
31,087; Clifford W. Browning, No. 32,201; R. Randall Frisk, No. 32,221; Daniel J. 
Lueders, No. 32,581; Michael D. Beck, No. 32,722, Kenneth A. Gandy, No. 33,386; 
Timothy N. Thomas, No. 35,714; Kerry P. Sisselman, No. 37,237; Kurt N. Jones, No. 
37,996; John H. Allie, No. 39,088; Jeffrey A. Michael, No. 37,394; Holiday W. Banta, No. 
40,311; Troy J. Cole, No. 35,102; L. Scott Paynter. No. 39,797; J. Andrew Lowes, No. 
40,706; Darrin Wesley Harris, No. 40,636; Matthew R. Schantz, No. 40,800; Gregory B. 
Coy, 40,967; Lisa A. Hiday, No. 40,036; John V. Daniluck, 40,581 and James L. Rowe, 
No. 18,448, all of the same address and who are members of the Bar of the State of 
Indiana, constituting the firm of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, Indianapolis, Indiana 46204 United 
States of America 

[X] as agent [ ] as common representative 

to act on his/her behalf before the competent International Authorities in connection with 
the international application concerning PARTICLE ANALYSIS SYSTEM AND METHOD, 
Agent's File Reference PUR-63/75P, International application Number unknown at this 
time, filed with the United States receiving office and to receive payments on his/her 
behalf. 

(Place) t l](TU6fcr\ ^ TX (Date) i^' [2^29^9^ 

Name: Jeff^ry KAO 
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PATENT COOPERATION TREATY 
Appointment of Agent or Common Representative 



The undersigned applicant, Huabei JIANG, hereby appoints Harold R. Woodard, No. 
16,214; C. David Emhardt. No. 18,483; Joseph A. Naughton, Jr., No. 19,814; John V. 
Moriarty. No. 26,207; John C. McNett, No. 25,533; Thomas Q. Henry, No. 28,309; James 
M. Durlacher, No. 28,840; Charles R. Reeves, No. 28,750; Vincent O. Wagner, No. 
29,596; Steve Zlatos, No. 30,123; Spiro Bereveskos, No. 30,821; William F. Bahret. No. 
31,087; Clifford W. Browning, No. 32,201; R. Randall Frisk, No. 32,221; Daniel J. 
Lueders, No. 32,581; Michael D. Beck. No. 32,722, Kenneth A. Gandy, No. 33,386; 
Timothy N. Thomas, No. 35,714; Kerry P. Sisselman, No. 37,237; Kurt N. Jones, No. 
37,996; John H. Allie, No. 39,088; Jeffrey A. Michael, No. 37,394; Holiday W. Banta. No. 
40,311; Troy J. Cole, No. 35,102; L. Scott Paynter, No. 39,797; J. Andrew Lowes, No. 
40,706; Darrin Wesley Harris, No. 40,636; Matthew R. Schantz, No. 40,800; Gregory B. 
Coy, 40,967; Lisa A. Hiday, No. 40,036; John V. Daniluck, 40,581 and James L. Rowe, 
No. 18,448, all of the same address and who are members of the Bar of the State of 
Indiana, constituting the firm of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, Indianapolis, Indiana 46204 United 
States of America 



to act on his/her behalf before the competent International Authorities in connection with 
the international application concerning PARTICLE ANALYSIS SYSTEM AND METHOD, 
Agent's File Reference PUR-63/75P, International application Number unknown at this 
time, filed with the United States receiving office and to receive payments on his/her 
behalf. 
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(54) Title: PARTICLE ANALYSIS SYSTEM AND METHOD 
(57) Abstract 

A system (20) and method are disclosed for the 
self-calibrating, on-line determination of size distribution f(x) 
and volume fraction <i> of a number of particles (P) dispersed 
in a medium (M) by detecting one or more propagation 
characteristics of multiply scattered light from the particles 
(P). The multiply scattered light is re-emitted in response to 
exposure to a light source (21) configured to provide light at 
selected wavelengths. The determination includes calculating 
the isotropic scattering and absorption coefficients for the 
particles (?) by comparing the incident and detected light to 
determine a measurement corresponding to the propagation 
time through the scattering medium (M), and iteratively 
estimating the size distribution f(x) and volume fraction $ 
as a function of the coefficients for each of the wavelengths. 
An estimation approach based on an expected form of the 
distribution and the mass of the particles is also disclosed. 
Furthermore, techniques to determine a particle structure factor 
indicative of particle-to-particle interactions which vary with 
particle concentration and influence light scattering at high 
concentrations is disclosed. 
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droplets dispersed in a fluid medium, the droplets and the medium having 
different indices of refraction. 
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(d) determining a calculated isotropic scattering coefficient for each 
of the wavelengths from an estimate of at least one of particle size distribution 
or particle volume fraction; 

(e) comparing the observed isotropic scattering coefficient and 
calculated isotropic scattering coefficient for each of the wavelengths to 
establish and error; 

(f) adjusting the estimate and repeating said determining and said 
comparing until the error reaches a desired minimum; and 

(g) providing an output corresponding to at least one of the particle 
size distribution, particle volume fraction, or a particle interaction parameter. 

47. The method of claim 46, wherein said determining includes 
establishing the particle interaction parameter as a function of the particle 
volume fraction and the particle size distribution. 

48. The method of claim 47, wherein the particle interaction 
parameter corresponds to the P-Y hard sphere structure factor. 

49. The method of claim 48, wherein the calculated isotropic 
scattering coefficient is determined with an equation having the structure 
factor and an estimated product of the size distribution and the volume 
fraction as arguments. 

50. The method of claim 46, wherein the particles include liquid 
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function of mass of the particles. 

43. The system of claim 41 , further comprising a second sensor 
providing a second detection signal, said scattering signals being determined 

5 as a function of said first and second detection signals for each of the 
wavelengths. 

44. The system of claim 41 , wherein said structure factor 
corresponds to a P-Y hard sphere structure factor model. 



45. The system of claim 41 , wherein said processor further 
generates an absorption signal corresponding to an absorption coefficient of 
the medium and determines said absorption signal from said propagation 
signal. 



46. A method of particle analysis, comprising: 

(a) exposing a number of particles to a number of light wavelengths 
of predetermined time-varying intensity; 

(b) detecting multiply scattered light from the particles in response 
20 to said exposing to determine a number of time-based propagation 

characteristics of the particles each corresponding to a different one of the 
wavelengths; and 

(c) calculating an observed isotropic scattering coefficient for each 
of the wavelengths from the characteristics; 



10 



15 
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different predetermined wavelengths of incident light each having a 
predetermined time-varying intensity; 

a first sensor spaced apart from said source, said first sensor being 
configured to provide a first detection signal corresponding to detected light, 
5 the detected light being multiply scattered by the particles; 

a processor responsive to said first detection signal and being 
configured to receive an exposure signal corresponding to said incident light, 
said processor being configured to generate: (a) a number of propagation 
signals by comparing said first detection signal and said exposure signal for 

io each of said wavelengths, said propagation signals each characterizing time 
of flight of the detected light through the medium resulting from multiple 
scattering by the particles for a corresponding one of said wavelengths, (b) a 
number of scattering signals each corresponding to an isotropic scattering 
coefficient of the medium and being determined from a corresponding one of 

15 said propagation signals, and (c) an output signal indicative of at least one of 
size distribution or volume fraction of the particles, said output signal being 
determined from said scattering signals and a structure factor, said structure 
factor accounting for particle-to-particle interactions influencing light scattering 
behavior of the particles; and 

20 an output device responsive to said output signal to provide an output 

corresponding to said size distribution or said volume fraction of the particles. 

42. The system of claim 41 , wherein said processor is further 
configured to determine said size distribution or said volume fraction as a 
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calculating a number corresponding to the isotropic scattering coefficient as a 
function of an estimate of the volume fraction and the size distribution. 

37. The method of claim 32, wherein the incident light is intensity 
5 modulated at a predetermined frequency, and said determining includes 
comparing the incident light and the scattered light to determine the value, 
and the value is representative of a relative phase or amplitude of the 
scattered light for the predetermined frequency. 

io 38. The method of claim 37, wherein said detecting includes 

detecting the scattered light with a second sensor spaced apart from the first 
sensor by a separation distance and said determining includes calculating the 
quantity in accordance with the separation distance. 

15 39. The method of claim 37, wherein said providing includes 

determining the volume fraction or size distribution in accordance with the 
diffusion equation for multiply scattered light. 



40. The method of claim 32, wherein the particle interaction 



20 



parameter is determined from a P-Y structure factor model. 



41 . A system for analyzing a number of particles suspended in a 



medium in sufficient concentration to multiply scatter light, comprising: 



a light source configured to expose the medium to a number of 



68 



WO 98/20323 




YUS97/20539 



10 



15 



(b) detecting multiply scattered light in response to said exposing to 
determine a time-based value characteristic of propagation time of the 
multiply scattered light through the fluid; 

(c) determining a quantity as a function of the value, the quantity 
corresponding to an isotropic scattering coefficient; and 

(d) providing an output determined from the quantity, the output 
corresponding to at least one of particle volume fraction, garticle size 
distribution, or a particle interaction parameter, the particle interaction 
parameter corresponding to particle-to-particie interactions influencing light 
scattering behavior of the particles. 

33. The method of claim 32, further comprising controlling a process 
in which the particles are altered by utilizing the output as a feedback 
variable. 

34. The method of claim 32, wherein said providing includes 
establishing an estimate corresponding to volume fraction and the size 
distribution and determining the particle interaction parameter as a function of 
the estimate. 

35. The method of claim 34, further comprising constraining the 
estimate to maintain mass balance. 

36. The method of claim 32, wherein said providing includes 
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corresponding to a nonlinear relationship between particle concentration and 
an isotropic scattering coefficient for the particles. 



27. The method of claim 26, wherein the particles are suspended in 
5 a fluid medium and further comprising controlling a process in accordance 

with the output. 

28. The method of claim 26, wherein the particles have a 
concentration in the fluid medium of at least about 10% by volume. 

10 

29. The method of claim 26, wherein said determining includes 
calculating at least one of volume fraction and size distribution of the particles 
in the fluid. t 

15 30. The method of claim 26, wherein the particle interaction 

parameter is determined from the Percus-Yevick hard_sphere model. 

31. The method of claim 30, wherein said calculating includes 
adjusting the hard sphere model to account for forces between the particles. 

10 

32. A method of analysis, comprising: 

(a) exposing a fluid to an incident light, the fluid having a number of 
suspended particles therein, the suspended particles being sufficiently 
concentrated in the fluid to scatter light; 
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error reaches a desired minimum. 

23. The method of claim 22, wherein said providing includes 
selecting the estimate to generally maintain mass balance of the particles, 
and the estimate corresponds to an expected form of the size distribution of 
the particles. 

24. The method of claim 22, wherein the calculated isotropic 
scattering coefficient is determined with an equation having an estimated 
product of the size distribution and Ihe volume fraction as an argument. 

25. The method of claim 19, wherein the particle interaction 
parameter is determined from a P-Y structure factor model. 

26. A method of particle analysis, comprising: 

(a) exposing a number of particles to a number of light wavelengths of 
predetermined time-varying intensity; 

(b) detecting multiply scattered light from the particles in response to 
said exposing to determine a number of values each corresponding to a 
different one of the wavelengths, the values each being representative of a 
time-of-flight characteristic of the particles; and 

(c) providing an output determined from the values, the output 
corresponding to at least one of a particle size distribution, particle volume 
fraction, or a particle interaction parameter, the particle interaction parameter 
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corresponding to at least one of a particle size distribution, particle volume 
fraction, or a particle interaction parameter, the particle interaction parameter 
corresponding to a nonlinear relationship between particle concentration and 
an isotropic scattering coefficient for the particles. 

20. The method of claim 1 9, further comprising controlling a process 
in accordance with the output. 

21. The method of claim 19, wherein said providing includes 
determining the observed isotropic scattering coefficient and an absorption 
coefficient for the particles for each of the wavelengths from a corresponding 
one of the values. 

22. The method of claim 19, wherein said providing includes: 

(i) determining an observed isotropic scattering coefficient for the 
particles for each of the wavelengths from a corresponding one of the values; 

(ii) establishing an estimate corresponding to at least one of the 
size distribution or the volume fraction; 

(iii) determining a calculated isotropic scattering coefficient for each 
of the wavelengths from the estimate; 

(iv) comparing the observed and calculated isotropic scattering 
coefficients to establish an error; 

(v) changing the estimate; and 

(vi) repeating said calculating, comparing, and changing until the 
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responsive to said output to regulate a reaction involving the particles. 

16. The system of claim 14, wherein said calculation means 
includes an estimating means for iteratively determining said size distribution 

5 or said volume fraction as a function of mass of the particles. 

17. The system of claim 14, wherein said structure factor is 
dependent on said particle size distribution and said particle volume fraction, 
and corresponds to a P-Y hard sphere model. 

io 

18. The system of claim 14, wherein said calculation means 
includes a number of parameters corresponding to a Weibull distribution. 

19. A method of particle analysis, comprising: 

15 (a) exposing a medium with a number of suspended particles to a 

number of light wavelengths, the wavelengths each being intensity-modulated 
at a predetermined frequency; 

(b) detecting multiply scattered light from the medium in response to 
said exposing to characterize propagation of the multiply scattered light 

20 through the medium with a number of values, the values each corresponding 
to a different one of the wavelengths and each being representative of at 
least one of a phase or an amplitude of the multiply scattered light relative to 
the predetermined frequency; and 

(c) providing an output determined from the values, the output 
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14. A system for analyzing a number of particles suspended in a 
medium in sufficient concentration to multiply scatter light, comprising: 

(a) a light source with a predetermined time-varying intensity 
configured to expose the medium to a number of different wavelengths of 
5 light; 



configured to provide a detected light signal corresponding to multiply 
scattered light from the particles at the different wavelengths in response to 
exposure to said source; 



including a calculation means for generating an output signal corresponding 
to at least one of particle size distribution, particle volume fraction, or a 
structure factor determined in accordance with an observed isotropic 
scattering coefficient of the medium, said structure factor being representative 

1 5 of particle interactions that influence light scattering behavior of said particles, 
said processor being configured to determine a value representative of said 
observed isotropic scattering coefficient from said detected light signal; and 

(d) an output device responsive to said output signal to provide an 
output corresponding to said at least one of said particle size distribution, said 

20 particle volume fraction, or said structure factor of the particles. 

15. The system of claim 14, further comprising a reaction vessel 
containing the particles and the medium, and wherein said output device 
includes a control element operatively coupled to said reaction vessel and 



(b) a sensor spaced apart from said source, said sensor being 



10 



(c) a processor responsive to said detected light signal, said processor 
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9. The system of claim 8, further comprising a reaction vessel 
containing the particles and medium, and wherein said output device 
includes a control element operatively coupled to said reaction vessel and 
responsive to said output to regulate a process involving the particles. 

10. The system of claim 8, wherein said processor is further 
configured to determine said size distribution or volume fraction as a function 
of mass of the particles. 



io 11. The system of claim 8, further comprising a second sensor 

providing a second detection signal, said comparison signals being 
determined as a function of said first and second detection signals for each of 
the wavelengths. 

15 12. The system of claim 8, wherein said processor further 

determines a structure factor indicative of particle-to-particle interactions, 
said structure factor varying in accordance with concentration of the particles 
in the medium. 



20 13. The system of claim 8, wherein said processor is configured to 

generate a particle interaction signal representative of particle-to-particle 
interactions that vary with particle concentration and influence light scattering 
behavior of the particles. 
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and a third value corresponding to mass of the particles. 

8. A system for analyzing a number of particles suspended in a 
medium in sufficient concentration to multiply scatter light comprising: 
5 a light source with a predetermined time-varying intensity configured to 

expose the medium to a number of different predetermined wavelengths of 
incident light; 

a first sensor spaced apart from said source, said first sensor being 
configured to provide a first detection signal corresponding to detected light, 
io the detected light being multiply scattered by the particles; 

a processor responsive to the first detection signal and being 
configured to receive an exposure signal corresponding to said incident light, 
said processor being configured to generate: (a) a number of comparison 
signals each corresponding to a difference with respect to time between the 
15 detected light and the incident light for a corresponding one of the 

wavelengths, (b) a number of scattering signals each correspondingly 
determined from the comparison signals and each corresponding to an 
observed isotropic scattering coefficient of the medium for a different one of 
the wavelengths, and (c) an output signal indicative of one of size distribution 
20 or volume fraction of the particles, said output signal being determined as a 
function of said scattering signals; and 

an output device responsive to said output signal to provide an output 
corresponding to the size distribution or volume fraction of the particles. 
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as a function of the estimate and determining the estimate to generally 
maintain mass balance. 



one of a number of parameters for an expected form of the size distribution. 

5. The method of claim 4, wherein the incident light includes a 
number of different wavelengths of light, the first value is determined for each 
of the wavelengths, the parameters correspond to a Weibull distribution, and 
said calculating includes: 

(i) establishing a number of particle size increments; 

(ii) determining a scattering efficiency and a mean cosine of 
scattering angle as a function of the wavelengths and the increments; and 

(iii) performing a first summation over a range of the particle sizes 
for each of the wavelengths, the first summation having the parameters, the 
scattering efficiency, and the mean cosine of scattering angle as arguments. 

6. The method of claim 5, wherein said comparing includes 
performing a second summation over a range of the wavelengths, the second 
summation including a numerical difference between the first and second 
values as an argument. 

7. The method of claim 5, wherein said changing includes updating 
the parameters as a function of the wavelengths, the first and second values, 



4. 



The method of claim 1 , wherein the estimate corresponds to 
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CLAIMS 



What is claimed is: 



A method of particle analysis, comprising: 
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(a) exposing a number of particles in suspension to an incident light 
with a predetermined time-varying intensity, the particles being sufficiently 
close to one another to multiply scatter light; 

(b) detecting multiply scattered light from the particles in response to 
the incident light to determine a first value corresponding to an observed 
isotropic scattering coefficient of the particles; 

(c) establishing an estimate corresponding to at least one of volume 
fraction or size distribution of the particles; 

(d) calculating a second value from the estimate, the second value 
corresponding to a calculated isotropic scattering coefficient; 

(e) comparing the first and second values to establish an error; 

(f) changing the estimate; and 

(g) repeating said calculating, comparing, and changing until the error 
reaches a desired minimum. 

2. The method of claim 1, wherein the incident light includes a 
number of different wavelengths of light, the first value is determined for each 
of the wavelengths, and said establishing, said calculating, said comparing" 
and said changing are performed for each of the different wavelengths. 

3. The method of claim 1 , further comprising controlling a process 
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embodiments have been shown and described and that all changes and 
modifications that come within the spirit of the invention are desired to be 
protected. 
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a fluid medium is determinable through measurements of multiplied scattered 
light in accordance with the present invention. Thus, the determination of 
liquid droplet size distribution to assess dispersion processes is a significant 
problem in industry that may be addressed by the present invention. 
5 In another embodiment, particle analysis techniques of the present 

invention include monitoring processes involving genetically engineered 
microbes. For example, the production of insulin in the granules of 
genetically engineered E. coli bacteria may be monitored to provide rapid 
feedback for optimization of the insulin production process. 
10 Also, it should be appreciated that the present invention finds 

application in many process involving the distribution of fine particles, 
including, but not limited to: processes involving a dispersion of liquid or solid 
particles in a fluid medium; chemical reactions having particles of mixed 
phases (such as solids and liquid droplets dispersed together in a common 
15 fluid medium); and the mechanical formation of particles. 

All publications, patents and patent applications cited in this 
specification are herein incorporated by reference as if each individual 
publication, patent, or patent application were specifically and individually 
indicated to be incorporated by reference. Copending U.S Patent Application 
20 Nos. 08/747,1 12, filed 8 November 1996, and 60/050,809, filed 26 June 1997 
are hereby incorporated by reference as if set-forth completely herein. While 
the invention has been illustrated and described in detail in the drawings and 
foregoing description, the same is to be considered as illustrative and not 
restrictive in character, it being understood that only the preferred 
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nanometers with the illustrated distributions being normalized to the peak 
value (vertical axis). Noise levels of 0%, 1%, and 5%, respectively, are 
shown in FIGS. 17a-17c. The results are tabulated in table 6 as follows: 





Actual 
Parameters 


0% Noise 


1 % Noise 


5% Noise 


Mean Diameter (nm) 


250 


250 


253.4 


243.9 


Percent Error (%) 




0.00 


1.36 


2.44 


Deviation (nm) 


50 I 50 

i 
I 


48.3 


32.9 


Percent Error (%) 




0.00 


3.4 


34.2 


Volume Fraction 


0.31 


0.31 


0.314 


0.291 


Percent Error (%) 




0.00 


0.127 


6.13 


Ar (nm) 


25 ! 25 


22.7 


28.6 


Percent Error (%) 


! 0.00 


9.2 


14.4 



5 



Table 6: Summary of Recovered Distribution for Soft Sphere Model 



This table indicates that the soft sphere model facilitates a good recovery of 
mean and volume fraction values despite the presence of noise. 
io Because the photon migration measurements of the present invention 

depend upon differences in refractive indices between the dispersed particles 
and the surrounding fluid, the size distribution of a liquid droplet dispersion in 
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straight, solid line corresponds to the linear relationship of isotropic scattering 
coefficients to particle concentration when particle interactions are ignored 
(the structure factor equal to unity, S(<7)=1); the circular symbols denote 
experimental measurements; and the dashed line indicates the nonlinear 
5 change of isotropic scattering coefficient with percent solids predicted by the 
P-Y hard sphere model. Notably, the P-Y model agrees well with the 
experimental data points, both illustrating a nonlinear relationship at higher 
concentrations. FIG. 16a shows the change in isotropic scattering coefficient 
with percent solids of the particles by volume at a wavelength of about 780 
io nanometers for a polydisperse polystyrene suspension with a mean particle 
diameter of about 0.4 micron. FIG. 16b shows the isotropic scattering 
coefficient at about 670 nanometers for a polydisperse polystyrene 
suspension having a mean particle diameter of about 0.15 micron. FIG. 16c 
illustrates the isotropic scattering coefficient at about a 670 nanometer 
1 5 wavelength for a polydisperse polystyrene suspension having a mean particle 
diameter of about 0.2 micron. It should be appreciated that the relationship 
between the isotropic scattering coefficient and particle concentration is . 
generally linear at low concentration levels, and becomes increasingly 
nonlinear as the concentration level increases. 
20 In yet another example, FIGS. 17a-17c illustrate the results of 

simulations of the soft sphere model for the structure factor. For these 
experiments, the inverse algorithm was applied to simulated measurements in 
the presence of various levels of noise. The particle size axis (the horizontal 
axis shown) was provided over a range between 15 nanometers and 1500 
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in accordance with the present invention for samples S1 and S2 
corresponding to FIGS, 12a and 12b, respectively. A log-normal Gaussian 
distribution was assumed for f(x) in both cases. The mean size and standard 
deviation for samples S1 and S2 are correspondingly presented in tables 4 
and 5 as follows: 





FDPM 


DLS 


X-ray 


Mean (|jm) 


0.310 


0.314 


0.318 


Std. Deviation (um) 


0.118 


0.112 


0.115 


Volume Fraction- 


0.0047 


0.0045 (By Evaporation) 





Table 4: Titanium Dioxide Suspension Sample S1 (See FIG. 12a) 





FDPM 


DLS 


X-ray 


Mean (urn) 


0.318 


0.375 


0.359 


Std. Deviation (urn) 


0.079 


0.143 


0.145 


Volume Fraction 


0.0044 


0.0045 (By Evaporation) 





10 

Table 5: Titanium Dioxide Suspension Sample S2 (See FIG. 12b) 



FIGS. 16a-16c comparatively illustrate the structure factor in terms of 
isotropic scattering coefficient (vertical axis) versus percent solids of particles 
15 by volume (horizontal axis). For each of FIGS. 16a-16c, the generally 
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Ti0 2 particles does not appear to appreciably detract from the relative 
accuracy of the PSD determination. 

FIGS. 12a-12b comparatively illustrate size distribution determinations 
using Frequency Domain Photon Migration (FDPM) analysis (solid line) of the 
5 present invention, Dynamic Light Scattering (DLS) (short dashed line), and X- 
ray scattering (long dashed line) for each of two samples of titanium dioxide 
suspensions, designated S1 and S2 t respectively. For the purpose of FDPM 
measurements, 1000 milliliters (ml) of each sample was placed in a 14.6 
centimeter (cm) cylinder glass vessel with an inside diameter of 1 1.1 cm. 
io Sinusoidaliy modulated light was deliver to, and collected from each sample 
via optical fibers of 1000 micrometers in diameter. Both the source and the 
detector optical fibers were placed parallel to the axis of the cylinder with their 
tips close to the center of the vessel. Phase shift and amplitude modulation 
was measured over a modulation frequency range of about 16 megahertz 
1 5 (MHz) to 240 MHz. Relative phase shifts were calculated corresponding to 
relative source-detector separations ranging form 0.1 cm to 0.5 cm. Samples 
were gently stirred between measurements to prevent any sedimentation. 
The solution of the diffusion equation (1) for the transport of photons in 
multiply scattering solution with infinite media boundary condition was applied 
20 to the data of relative phase shifts at different modulation frequencies. Upon 
applying Marquardt-Levenberg nonlinear regression analysis to this data, 
isotropic scattering coefficients (p s ') were obtained at seventeen different 
wavelengths ranging from 400 nanometers (nm) to 900 nm. Particle size 
distribution, f(x), and volume fraction, <J> , were determined from these values 
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from experimental measurements for three different suspensions with the 
additional constraint of mass balance. The weighting parameter used for 
obtaining these results was a constant of 0.1. It has been discovered that a 
weighting parameter in the range between 0.01 and 0.2 works well for these 
samples. Table 3 presents the results from the volume fraction calculations 
for the three different suspensions as follows: 





Values of from Photon 
Migration (%) 


Values of 4> from DLS (%) 


PP722 


0.98 


0.95 


PP755 


1.61 


1.53 


PP788 


0.67 


0.63 



Table 3 Solids volume fractions for suspensions PP722, PP755 and 
io PP788 obtained with mass balance constraint. 



FIG. 12 illustrates a Photon Migration (PM) measurement of multiply 
scattered light to experimentally determine PSD (solid line) for an aqueous 
Ti0 2 suspension following a procedure similar to the polystyrene suspensions 
1 5 described above. For comparison, size distributions obtained from Dynamic 
Light Scattering (DLS) (dashed line) and determined from Transmission 
Electron Microscopy (TEM) (circle symbols) are shown. It should be noted 
that the size distributions obtained from photon migration measurements 
straddled those obtained from DLS and TEM. The non-spherical nature of 
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for each inverse solution required a few seconds on a SunSparc 10 
Workstation using the data gathered by the IBM compatible system. The 486 
based IBM compatible computer and SunSparc 10 Workstation collectively 
provide one example a processor suitable for use with the present invention; 
5 however, in other embodiments, a single piece of dedicated equipment, or 
other arrangement is envisioned as would occur to one skilled in the art. 

FIGS. 6 through 8 show reconstructed particle size distributions from 
experimental measurements for the three different suspensions without the 
mass balance constraint. The solid lines (PM) denote the reconstructed 
io values from the photon migration measurements of multiply scattered light 
while the dashed lines represent the size distribution obtained from DLS for 
comparison. Table 2 presents the volume fraction determination for the three 
different suspensions, as follows: 





Values of <|> from Photon 
Migration (%) 


Values of <j> from DLS (%) 


PP722 


1.16 


0.95 


PP755 


1.82 


1.53 


PP788 


0.75 


0.63 



15 

Table 2 Solids volume fractions for suspensions PP722, PP755 and 
PP788 obtained without mass balance constraint. 

FIGS. 9 through 1 1 display the reconstructed particle size distributions 
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Figure 4 illustrates the typical measurements of relative phase-shift, 
6 re |, on 0.96%. by volume PP722 sample as a function of modulation 
frequency at a wavelength of 800 nm. Measurements were conducted with 
source-detector separations |r $ -r d | of 2.0, 1.5, and 1.0 cm. Consequently, 
5 relative phase-shift was reported at relative detector-detector separations |r d1 - 
r d2 | of 1 .0 and 0.5 cm. It is noteworthy that the two sets of phase-shift data 
recorded for each of the 0.5 cm separations were approximately the same, 
and the relative phase-shift for the 1 .0 cm separation was about twice as 
large as the relative phase shift corresponding to the 0.5 cm separations. 
io This dependence of phase shift on the relative separation of the source and 
detection sites is consistent with presence of the |r d1 -r d2 | term in equation (7). 

FIG. 5 graphically depicts experimental values (x symbols) for isotropic 
scattering coefficients obtained for the PP722 sample from measurements of 
multiply scattered light at 10 wavelengths between 360 and 800 nm. 
1 5 Expected values (solid line) predicted by Mie theory were calculated with 
equation (3) from experimental volume fraction and PSD. These 
measurements were conducted without the use of any external reflectance or 
calibration standard. 

The f(x) and <J> for each of the three polydisperse suspensions were 
20 determined assuming a Weibull distribution both with and without the mass 
conservation constraint. The initial estimates of the reconstruction 
parameters a,b t c,d were up to 100% greater than the output parameters. The 
particle size distributions, f(x), and volume fraction, $ , required 5 iterations 

until the function X" reached a desired minimum value. Computational time 
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detected with the reference PMT. In order to obtain values of relative phase- 
shift, 6rel» measurements were conducted as the detecting fiber optic was 
moved with micrometer precision to distances of 1.0, 1.5, and 2.0 cm away 
from the source fiber optic. At each source-detector separation, 
5 measurements were conducted at 10 modulation frequencies ranging from 24 
to 240 MHz. Relative phase-shift between any two source-detector 
separations was determined by subtracting the phase-sfiift measured at the 
farthest separation from that measured as a closer separation. By fitting 
experimentally measured relative phase-shift at varying modulation 
io frequencies to equation (7) using Levenberg-Maraquardt non-linear least 
square regression, parameter estimates of pa and p s ' were obtained. 

Three polydisperse samples of polystyrene microsphere suspensions 
were characterized using frequency domain measurements of multiply 
scattered light with this instrumentation. The three polystyrene samples were 
1 5 obtained from DOW Chemical Company of Midland, Michigan and designated 
PP722, PP755, and PP788 with mean particle diameters of 0.5763, 0.2954, 
and 0.1685 microns, respectively, as determined by independent DLS 
measurements (Microtrac Ultrafine Particle Analyzer, Honeywell, Leeds and 
Northrup, St. Petersburg, FL). The concentration of solids for PP722, PP755, 
20 and PP788 were 47.95%, 52.74%, and 46.03% solids by volume respectively 
as determined through evaporation measurements. The samples were 
diluted about 50 times with deionized ultrafiltered water to obtain samples of 
approximately 1% solids by volume. This concentration results in an opaque 
solution suitable to multiply scatter light. 
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outer diameter and 8 inch height. The tank was filled with a sample having 
particles suspended in a fluid medium. 

The light which propagated from the source fiber within the multiple 
scattering medium was collected by another fiber held vertically in place in the 
5 tank. The detected light was then delivered to a sample sensor, also a PMT, 
via an optic fiber. Fourier analysis of the 4 MHz pulse train delivered to the 
sample and reference PMT's yielded a series of harmonic intensity- 
modulation frequencies at multiples of 4 MHz, which upon application of 
heterodyne detection techniques permitted isolation of the individual intensity- 
io modulated frequencies. Heterodyning was performed by gain modulating the 
reference and sample PMTs at a harmonic of the laser repetition rate plus an 
offset cross-correlation frequency of 100 Hz. Gain modulation was 
accomplished using a commercial electronics package (ISS, Champaign, IL) 
and a -3 dBm RF signal from a frequency synthesizer (Marconi Instruments 
1 5 Signal Generator 2022A) amplified by a power amp (model 1403LA, ENI, 

Rochester NY). The resulting 100 Hz electronic signaLfrom the heterodyned 
PMTs provided the frequency domain phase-shift and amplitude modulation 
information. A conventional data acquisition module (ISS, Champaign, IL) 
was employed in a 486 IBM compatible personal computer to collect the 
20 measurement data. 

By sweeping the RF signals at varying harmonics of the laser repetition 
rates, measurements of phase-shift 9(co) and amplitude demodulation M(co) of 
a multiply scattered photon density wave propagating at varying modulation 
frequencies were detected with the sample PMT relative to the incident light 
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EXPERIMENTAL SECTION 



The present invention will be further described with reference to the 
following specific examples. It will be understood that these examples are 
illustrative and not restrictive in nature. 

Experimental measurements to determine PSD and volume fraction 
where conducted with a instrumentation configuration substantially similar to 
system 20. This instrumentation included a 10W argon ton laser (Beam-Lok 
2060, Spectra Physics, CA) which pumped a picosecond pulsed 
Titanium:Sapphire laser (model 3950B Tsunami, Spectra Physics, Mountain 
View, CA). The Ti:Sapphire laser produced an optical pulse train of equally 
spaced light pulses of 2 psec FWHM at a repetition rate of 80 MHz. The 
output beam was sent to a pulse picker (model 3980, Spectra Physics) having 
an AOM crystal to produce a pulse repetition rate of 4 MHz at 2 ps FWHM 
and a SHG crystal to provide for light wave frequency doubling (wavelength 
halving). As a result, wavelengths from 720 to 900 nm, and from 360 to 450 
nm were available for selection. The average power exiting the AOM crystal 
and the frequency-doubler (SHG crystal) was 50 mW and 10 mW, 
respectively. Depending upon the wavelength desired, the pulsed laser beam 
from the AOM or frequency doubler was sent to a glass slide beam splitter 28 
that directed approximately 20% of the light to a reference sensor, 
photomultiplier tube (PMT) (Hamamatsu R928, Hamamatsu, Japan), via a 
1000 pm optical fiber (HCP-M1000T-08, Spectran, Avon, CT). The remaining 
light was directed to the scattering sample by another optical fiber which was 
held stationary and vertical in a clear acrylic cylindrical tank having a 9 inch 
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radius that accounts for the degree of particle interactions. Moreover, 
besides structure factor, other particle interaction parameters may be 
selected and applied in accordance with the present invention to account for 
variation in light scattering properties with particle concentration due to 
5 particle-to-particle interactions. Such particle interaction parameters generally 
characterize nonlinearity between a scattering property of the particles, such 
as the isotropic scattering coefficient, and the particle concentration level. 
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art. 

In one such alternative structure factor model, the hard sphere model 
is modified to account for electronic affects beyond the particle surface. This 
phenomenon is particularly substantial for lower concentrations of 
5 electrolytes, where particle interactions are still significant. FIG. 14 

comparatively represents this modification as a soft sphere SS next to hard 
sphere HS. Soft sphere SS has a "soft" outer shell region corresponding to 
the interactive reach of the particle beyond its surface. This "soft" region. is 
designated by the radial extension increment Ar and accounts for the 
io electronic forces extending beyond the surface of the particles. The effective 
electrostatic charge EC2 is shown as a broken line at r + Ar. The sum r + Ar 
is treated as the "equivalent" hard sphere radius for the P-Y model which is 
otherwise applied in the same way as described above. FIG. 15 depicts the 
relationship of the isotropic scattering coefficient to light wavelength for the 
1 5 soft sphere model (diamond-shaped data points), the hard sphere model 
(triangular-shaped data points), and the unity structure factor, S(qr) = 1 
(square-shaped data points). 

Regardless of the structure factor model selected, any of the 
embodiments previously described and any such variations as would occur to 
20 those skilled in the art may be adapted to incorporate the structure factor 
technique. By way of nonlimiting example, the mass balance constraint 
described in connection with equations (22)-(30) or the Weibull distribution 
simplification described in connection with equations (15)-(20) may be 
included in a process to determine a structure factor or effective hard sphere 
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which isotropic scattering coefficients for each wavelength (i^ s »)j are 

calculated through the summation approximation of equation (31) for all 
wavelengths using the estimates of f(x) and <j>, and the corresponding 
calculated values of S(g), g, and q^. After stage 642, control flows to stage 
5 644 to determine x 2 a s described in connection with stage 544 of process 
500. The x 2 of observed and calculated coefficients is then compared to an 
error threshold in conditional 646 to determine if the current estimates of f(x) 
and 4> are suitable for output. If conditional 646 is satisfied, control flows to 
stage 650 to output f(x), <j>, and S(q). If conditional 646 is not satisfied, 
io control flows to stage 648. 

In stage 648, the estimation of f(x) and $ is updated. This update also 
includes recalculation of g and q scat based on the current S(q) value. Once 
updates are performed in 648, loop 640 closes back on stage 642 to provide 
updated calculations of the corresponding isotropic scattering coefficients. 
15 Loop 640 is iteratively executed until the x 2 difference is minimized to an 

acceptable level. The particle size distribution f(x), the volume fraction <j>, and 
structure factor S(q) may be output in stage 650. Aside from the improved 
accuracy of size distribution estimation and volume fraction at high particle 
concentrations, the structure factor S(q) is useful in ascertaining information 
20 about the behavior of colloidal suspensions and emulsions. For example, the 
structure factor alone is useful in monitoring colloidal suspensions and 
predicting shelf-life of emulsions with respect to phase separation. In other 
embodiments, different particle structure factor models may be employed 
other than the P-Y hard sphere model as would occur to those skilled in the 
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1 =1+ 24^ 



A(sin x - xcos x) + B {^~~ ' j 
^f24 ( 6 ) . r, 12 24 ^ 



x cos .t + 2 sin -x 

x 



COS X 



(32) 



where: a is the diameter of the particle determined from f(x); x=qra; 
cp=(7rNDa 3 /6); ND=number density of particles determined from a and <|>; 
5 A=[(1+2<p) 2 /(1-cp) 4 ]; and B=[-3cp(cp+2) 2 /2(1-<p) 4 ]. It should be noted that S(q) is 
a function of f(x) and <j>. 

FIG. 13 is a flow chart of process 600 which may be used to determine 
f(x) t <p ( and S(q). Process 600 may be performed using system 20, system 
120, or such other system as would occur to one skilled in the art. Stages 
io 510-520 of process 600 are the same as the like numbered stages of process 
500 described in connection with FIG. 2. After stage 520, stage 630 of 
process 600 is encountered. In stage 630, f(x) and <f> are estimated and S(q) 
is calculated. Notably, the determination of S(g) is dependent on f(x) and $ 
estimates for this embodiment. The S(q) calculation utilizes the mean particle 
!5 diameter for a as determined from the f(x) estimate, and the number density 
ND as determined from the f(x) and volume fraction <|> estimates. Also, g and 
Qscat are initially determined to fit the calculated S(q) to the measured isotropic 
coefficient fi s ' and absorption coefficient ji a in accordance with a summation 
approximation of equation (31) for each wavelength X r 
20 Next, loop 640 is entered. In loop 640, stage 642 is encountered in 
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A. = ITT J 2;r ^. ^ r 2 0 - r))S(q)<f>{\ - g(n, X, x))sm{e)d0 



(31) 

where: <7 = equals [^Trn/A^sin^))] and q scat is the angular scattering 
5 efficiency. The structure factor S(q) accounts for nonrandom particle 

structure within the fluid that typically varies with concentration of the particles 
in the medium. For example, the electrostatic charge on suspended particles 
generally influences spacing and organization of substantially concentrated 
particles in a fluid, which in turn influences light scattering. For low 
io concentrations, these forces are usually negligible, so the particles may be 
treated as random scatterers with S(q) = 1. For a structure factor of unity 
(S(q) = 1), the calculation of isotropic scattering coefficients m s ' at each 
wavelength of interest is the same in equations (3) and (31). For higher 
concentrations, S(g) typically varies with the structure imparted by the particle 
1 5 interactions. 

One approach to ascertaining an appropriate S(q) is to treat particles 
as hard spheres of radius r as represented in FIG. 14 by particle hard sphere 
HS. For a hard sphere, the electrostatic charge resides on the surface, as 
represented in FIG. 14 by broken line EC1. For this approximation, the 
20 Percus-Yevick (P-Y) hard sphere model may be employed to predict S(g) as 
follows: 
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boundary conditions. In other embodiments, the present invention may be 
readily adapted to different boundary conditions. For example, in one 
alternative embodiment, a processing chamber or passageway containing the 
scattering media may be interrogated by detecting multiply scattered light that 
5 is reflected or backscattered by the media through a "viewing window" 
arrangement instead of directly immersing fibers in the media. This 
arrangement results in a reflectance boundary condition. The present 
invention may be adapted using techniques known to those skilled in the art 
to accommodate this reflectance geometry or such other boundary conditions 
10 as may arise from a particular probe arrangement. 

When the concentration of particles suspended in a medium is low 
(e.g. <2%), the particles act independently with regard to light scattering; 
however, as the concentration increases, interactions between suspended 
particles become more significant. These particle-to-particle interactions may 
15 impact accuracy of the f(x) and 4> determination in accordance with the 

present invention. Thus, a particle interaction parameter which accounts for 
concentration-dependent particle interactions that influence light scattering 
properties is desirable. Preferably, the determination of this particle 
interaction parameter may be incorporated into techniques of the present 
20 invention for determining f(x) and <(>. 

For one preferred embodiment, a structure factor S(q) for the particles 
has been found to be suitable as the particle interaction parameter. In this 
embodiment, the structure factor S(q) is introduced by modifying the isotropic 
scattering coefficient (m 3 ') determination of equation (3) as follows: 
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intensity. In other embodiments of the present invention, the light source 
arrangement is configured to provide wavelengths and time-varying intensity 
suitable for interrogation of the particular process being monitored using 
techniques known to those skilled in the art. 
5 Analyzer 150 also includes optical fiber 134 to detect re-emitted light 

from the particles in vessel 160 after multiple scattering. It is preferred that 
analyzer 150 be of rugged design suitable for industrial applications. 
Analyzer 150 includes a processor configured to execute software to 
determine values representative of PSD or volume fraction in accordance 
io with the present invention. Furthermore, analyzer 150 is also operatively 
coupled to control elements 170 which are used to regulate processing in 
vessel 160. By way of non-limiting example, elements 170 may include 
valves, heaters, or agitation devices electronically regulated by analyzer 150. 
Analyzer 150 is further configured with appropriate programming and 
1 5 interfaces to provide one or more output signals to control elements 170 as a 
function of PSD or volume fraction. Thus, system 120 provides a closed loop 
feedback control system capability for regulating chemical processes from on- 
line measurement of PSD and/or volume fraction. 

Systems 20 and 120 of FIGS. 1 and 3, respectively, are described with 
20 corresponding light source fibers 30, 130 and detecting fibers 36, 136 

geometrically arranged to approximate infinite boundary conditions — such 
that light travelling from fibers 30, 130 to fibers 36, 136 does not encounter a 
substantial boundary. Equations (6) or (7) are applied to determine the 
isotropic scattering coefficient p s ' and the absorption p a under such infinite 
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This mass balance constraint serves as an additional regularization scheme 
5 and typically provides further stabilization. 

FIG. 3 depicts system 120 of an alternative embodiment of the present 
invention is depicted. System 120 includes analyzer 150 coupled to reaction 
vessel 160 by fibers 130 and 134. Analyzer 150 is configured to determine 
PSD and volume fraction for particles dispersed in a fluid medium in vessel 
10 160 in accordance with the present invention. Generally, the optimum 

number and size of the wavelengths selected for interrogation depends on 
the nature of the particles, medium, and process involved. Furthermore, it 
has been discovered that typically fewer wavelengths may be required to 
determine volume fraction as compared to PSD for a given arrangement. 
15 Analyzer 150 includes a light source configured to deliver light to 

vessel 160 via fiber 130. This light source may resemble source 21 of system 
20, but preferably is a lamp-based or laser diode array system configured to 
provide light of a desired time-varying intensity at selected wavelengths. In 
one embodiment, system 120 includes an array of laser diodes which each 
20 correspond to a desired wavelength suitable for the particular process and 
material being interrogated in vessel 160. For this embodiment, the diodes 
are configured to provide incident (output) light with a sinusoidally modulated 
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x 2 = Z[(m;);-(k);] 2 + w[m, - r P — *xV(x)dxT 

L ■» 12 J (23) 

where w is a weighting parameter. To efficiently update the estimate 
and account for the relationship of equation (23), the following Jacobian- 
based relationship is employed: 

(3 r 3 + + <xI)Ac = 3 T [(,a $ y - (^) c ] + V (24) 

where 
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,5 " L ' •» H 12 * 1* y l2 dc (29) 
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particle sizes which surround the pth particle size. Current experiments have 
shown c=0.35 to be one preferred condition for the filter. This low pass filter 
typically enhances the particle size distribution recovery considerably by 
relaxing the requirement of the selected type of size distribution. It should be 
appreciated that the wider the distribution, the larger the required filter width 
for enhanced reconstruction. From the "filtered" f(x), the calculated isotropic 

scattering coefficient (m s -)^ is computed for updating the Jacobian matrix. 

In still a further embodiment, calculation of the inverse solution is 
improved when the total mass of the particles M t and the mass density p are 
known. For this embodiment, the estimation of f(x) and <)> focuses not only 
on the comparison of observed and calculated isotropic scattering 
coefficients, but also compares the known total particle mass to the mass 
obtained as a function of p and the estimates as follows: 

M t = /P-Ux 3 [<t>f(x)]dx. (22) 
o 12 

15 For some process configurations "mass balance" may often serve as a 

process control parameter in addition to PSD and volume fraction. For 
example, in emulsion polymerization processes, the total mass of monomer 
dispersed prior to initiation of polymerization is known. Proper reconstruction 
of the size distribution of the dispersed monomer phase should enable 

20 prediction of the total mass of dispersed monomer. Disappearance of the 

mass of monomer should reappear as mass of polymer. To include this mass 
balance constraint in the algorithm, we consider it as a constraint in the least 
squares minimization comparison as follows: 
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L 3a 3b 3c 3d (18) 

The elements of the Jacobian are numerically computed using 
equation (3). The vectors, (Ms') 0 and (Ms') c contain the observed and 
computed values of isotropic scattering coefficients at the wavelengths of 
5 interest: 

(19) 

A£ is the vector updating the four parameters, a, b, c and d: 
A<; = (Sa,Sb,5c,Sd) T (20) 
io Thus, the particle sizing task now becomes to recover the four parameters (a, 
b ( c and d) to describe f(x) and 4> using the approach described in connection 
with loop 540. 

In still another embodiment, a low pass filter is employed to smooth the 
estimate based on the parameters for an expected form of distribution, such 
15 as the Weibull distribution. In this embodiment, the size distribution is 

computed from the parameters and then subjected to a digital low pass filter 
by averaging over a window of N*Ax. 

f(x p )- = (i -q)f( Xp r P r^.)" IJ 

,=P "T (21) 
where 4 is a factor between 0 and 1 , and the summation is over the N* 



WO 98/20323 



PCT/US97/20539 



10 



assuming an expected form of the PSD. Generally, particulate processes can 
usually be characterized by Gaussian or log-normal distributions, both of 
which can be represented by proper parametric choices for a Weibull 
distribution described as follows: 

J fx-aV"l 

(15) 

where a, b and c describe the peak location, width, and shape, respectively, 
of the distribution. For this embodiment, d corresponds to <J> , and the 
equation (3) summation approximation becomes: 



( bl_ b~~J Jj 



^i!^^'-^ ))J 



cl C - 
b 



r 



-a 



exp 



(16) 



where the expected range of particle diameter x from x^to x^ is determined 
and incremented into k number of discrete sizes x t = x m(n + (i-1) Ax, 
Ax = (xma* - xmin) / (z - 1) for i=1 to z, and the terms QscaAj.n.X;) and g(\ f n t x t ) 
1 5 are determined from look-up tables based on a^, n, x r 

Using Newton's method, the updates based estimation of the a,b,c,d 
parameters can be obtained from the resulting system of equations: 

3 T 3A<;=3 T [( M :r-ai;n 

where the Jacobian matrix 3 now represents the sensitivity of isotropic 
20 scattering coefficients measured at wavelengths j=1 ,M upon the four 
parameters (a, b, c and d): 
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photon quantity in terms of intensity l(t), and the horizontal axis represents 
time t. The TOF distribution may be employed to obtain the observed 
isotropic scattering coefficient by applying the time domain form of the 
diffusion equation. Typically, time of flight values are in the range of 
5 picoseconds to a few nanoseconds. 

In comparison, the measurement of intensity modulated light 
propagation in the frequency domain in accordance withrequations (1)-(7) and 
accompanying text is provided in terms of phase shift 9 and amplitude 
modulation M in FIG. 2b. Phase shift d is represented by the left vertical axis, 

10 amplitude modulation M is represented by the right vertical axis, and 

frequency is represented by the horizontal axis which increases from left to 
right. The relationship of frequency to phase shift is characterized by the 
continuous curve that increases with increasing frequency (as designated by 
the left-pointing arrow) and the relationship of frequency to amplitude 

15 modulation is characterized by the continuous curve which decreases with 
increasing frequency (as designated by the right-pointing arrow). 

It should be appreciated the "time-of-flight" characteristics of a 
scattering medium alternatively may be described in terms of phase shift and 
amplitude modulation relative to the modulation frequency. Notably, in both 

20 the frequency domain and time domain approaches, an intrinsic optical 

characteristic or optical property of the particles, such as isotropic scattering 
and absorption coefficients, may be determined from a time-based 
propagation measurement of multiply scattered light relative to input light. 

In another embodiment, the determination of f(x) and 4 is simplified by 
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Once the estimate is updated in stage 548, control returns to stage 
542, via loop 540 to calculate the isotropic scattering coefficients based on 
the updated estimate. Repeated executions of loop 540 are anticipated 

before the desired minimum value of x 2 is obtained. 
5 Preferably, stages and conditionals corresponding to reference 

numbers 520 to 550 are performed through a software routine executed by 
processor 50, and processor 50 is of a type suitable to perform the extensive 
calculations of process 500 in an adequate amount of time. Notably, the 
phase shift may be determined from detection at a single location rd using 

10 equation (6) or at multiple locations to determine the 0rel as a function of rd1 
and rd2 using equation (7). Furthermore, instead of phase shift, process 500 
may be adapted to determine f(x) and <{> based on amplitude M 
measurements in the frequency domain. 

Referring back to FIG. 1, it should be appreciated that the length of 

is scattering path SP is longer than the straight, direct pathway from fiber 30 to 
fiber 36. This pathlength increase due to the scattering media also results in 
a comparatively longer propagation time or "time of flight." Because photons 
typically travel along many different paths in a scattering media, the varying 
times of flight of the photons is usually amenable to description in a time of 
20 flight distribution. Thus, in one alternative embodiment, the diffusion equation 
may be utilized to compare incident light to resulting multiply scattered light in 
the time domain. For example, FIG. 2a illustrates the relationship between an 
input light pulse LP of a selected wavelength to a corresponding "time of 
flight" response, shown as distribution TOF; where the vertical axis represents 
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3 T 3AC=3 T [(^;r-(^) c ] 



(11) 

which results in a system of linear equations that may be solved using 
Newton's method to obtain the vector AC, containing the update values of 
(|>f(x) across all bins, h=1, N as follows: 

5 AC = (S((()f(x)) l .5((})f(x)) 2 ,...5(<(>f(x)) N ) T (12) 

The vectors, (ja s )° and (ns) c contain the observed and computed 
values of isotropic scattering coefficients at the wavelengths of interest, j=1,M 
are: 

(M;)° = ((rt)^(K)r,.»-(rt)L) ( (K) e =((^)^(rt)^.-.(n;)^) (13) 

io The matrix 3 T 3 in equation (11) is sometimes ill-conditioned and may 

be more readily solved through regularization of the decomposition of 3 T 3. 
One type of regularization scheme is the Maraquardt-type which may be 
stated as: 

(3 T 3 + aI)A^=3 T [(^;)"-(^) c ] 



(14) 



where I is the identity matrix and a may be a scalar or a diagonal matrix. By 
adding a contribution to the diagonal terms in equation (14), we make 3 T 3 
more diagonally dominant which improves its invertability. The parameter a is 
usually determined by trial and error. Jiang, H. ( Paulsen, K., Osterberg, U., 
20 Pogue, B. and M. Patterson, Optical image reconstruction using frequency- 
domain data: simulations and experiments , J. Opt.Soc. Am. A 13, 253-266 
(1996) is cited as a source of additional information concerning these 
calculations. 
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where it is desired to minimize x so that it is less than or equal to an 
established error amount. Conditional 546 tests whether this minimization 
has taken place. If the minimum has been reached, then control flows to 
stage 550 and the estimate is provided as the measured product 4»f(x). 
Subsequently, f(x) and <J> may be separated by performing a summation 



However, if conditional 546 is not satisfied, then control flows to stage 
548 to improve the estimate of <J>f(x). One embodiment for rapidly converging 
to the desired minimum error is through application of a Jacobian matrix 
relationship. Specifically, the product <j>f(x) is divided into N number of bins 
representing <j>f(x) h . where h = 1 to N. In order to update <|>f(x) h in each bin, 
the Jacobian matrix which describes the sensitivity of the isotropic scattering 
coefficient to changes in <t>f(x) h in each bin h and at each wavelength ^is 
computed. This Jacobian matrix is given by: 



and each element is computed through a summation approximation of 
equation (3) with <frf(x) h and (|>f(x) h + A<j>f(x) h . The updates, 5((>f(x) h , 
resulting from the differences between measured and computed isotropic 
scattering coefficients, can then be obtained from: 



approximation of J<)>f(x)dx = <f> . 



o 




(10) 
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Once (m s ')^ has been obtained in stage 520, stage 530 is 
encountered. Stage 530 initiates an algorithm to inversely solve equation (3) 
for the integrands f(x) and <|> given(|a s «)° . This inverse solution starts by 
estimating the product <|>f(x)in stage 530. Once this product is known, f(x) 

TO 

5 and <|> may be separated by the relation J<f>f(x)dx = <j> . 

o 

Next, iteration loop 540 is entered at stage 542. In stage 542, the 
isotropic scattering coefficients for each of the wavelengths of interest are 
calculated from equation (3) using the estimate of 4>f(x). For this calculation, 
equation (3) is approximated by the summation equation (8) as follows: 

10 

_ f _ g(Xj .„, Xi)) Wxto . (8) 

3 i = l2 Xt X 

For equation (8), an expected range of particle diameter x from x min to 
x max is determined and incremented into k number of discrete sizes x t = x mjn + 
(i-1) Ax, where Ax = ( Xm ax - x m ,„) / (z - 1) for i=1 to z; and the terms 
15 Q*at(Aj. n . x i) and g^.n.xj are determined from look-up tables based on A 1 , n t 

Xj. The resulting set of calculated isotropic scattering coefficients, (^ s -)j are 

then compared to the observed isotropic scattering coefficients, (n s ')° in 
stage 544. This comparison is performed in accordance with: 

z 2 = Zk^);-(m;)?] 2 

(9) 
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Equations (6) and (7) may be expressed in terms of the absorption 
coefficient ua and isotropic scattering coefficient ns'=( 1 -g)ns through the 
substitution of equation (2) for D(/.). Therefore, given 6(rd) or 9 r el =|9(rdl)- 
9(rd2)l from the execution of loop 510, equations (6) and (7) may be fit, 
5 respectively, for an observed absorption coefficient an d isotropic scattering 
coefficient us' by applying a Levenberg-Maraquardt non-linear least square 
regression of the type disclosed in Press et al. f Numerical Recipes: The Art 
of Scientific Computing (Cambridge University Press 1992) to the data 
tabulated, as a result of the execution of loop 510 for each of the wavelengths 
io of interest. Preferably, processor 50 accumulates the frequency domain 
information in memory and is programmed to execute a Levenberg- 
Maraquardt non-linear least square regression algorithm to provide these 
optical coefficients. The resulting observed isotropic scattering coefficients 

for the wavelengths \ are represented as (ja s .)" . Preferably, (1^')° is 

\5 obtained by applying equation (7) to G r el measurements in stage 520. 

Unlikexonventional PSD measurement systems, these frequency 
domain measurements offer the ability to monitor absorption and scattering 
properties in multiply scattering systems without the need for spectral 
calibration of lamps or detectors to an external standard. Also, because the 

20 absorption and isotropic scattering coefficients do not appear as a product in 
equation (7), regression of phase-shift data provides an efficient means of 
determining optical characteristics from measurements of multiply scattered 
light. 
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10 



15 



position, r<j: 



9(r d ,co,A.) = tan 



„ 1 ImO( r tl ,co,A.) 
Re O(r d ,03 r A) 



(4) 



and 



iV/(r d ,<y,A) = ^j[\m<X>(r d . co, /I)] 2 + [Re *(r d , a;, /I)] (g) 

The solution to equation (1) can be obtained for reflectance or 
transillumination measurements on semi-infinite random media as well as for 
finite geometry's in which the source and detector are separated by distance 
|r s -rd|. For simplicity, an infinite geometry is assumed which results in the 
following approximation: 



(3///U){(/i,U)^) :+ ^ 3 } T 



sin 



tan 



CO 



(6) 



Measurement of phase-shift usually includes a contribution due to the 
electronics and photodetection systems. Upon conducting a spatially resolved 
phase-shift measurement at two or more consecutive detector positions, rd1, 
rd2, instrument errors are accounted for. The relative phase, Opel can then 
be predicted from the following relationship: 



CO' 



sin| — tan 



(O 



(7) 
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coefficients, s, and concentration, [C] f for all light absorbing constituents 
present; and D(X) is an "optical diffusion coefficient" which is given by: 

Da> = ! 

3|>i a (>o + n -gJhW]. (2) 

The following are cited as additional sources of information concerning 
5 these relationships. A. Ishimaru, Wave Propagation and Scattering in 
Random Media , Academic Press, New York (1976); S. Chandrasekhar, 
Radiative Transfer , Oxford University Press, New York (1960); and R. C. 
Haskell, Svaasand, L.O., Tsay, T-T, Feng, T-C, McAdams, M.S., and B.J. 
Tromberg, Boundary Conditions for the Diffusion Equation in Radiative 
io Transfer , J. Opt. Soc. Am. A 11, 2727-2741 (1994). 

The term (1-g)Ms = Ms' is the isotropic scattering coefficient which 
arises from multiple scattering and is dependent upon the particle size 
distribution, f(x) (where x is the diameter of the particles), and the total volume 
fraction, <j> as follows: 

^Qenr (x.n.A.)fI - u(x.n.A,)l ir/ 
15 o ~ x (3) 

where g is the mean cosine of the scattering angle from a single particle and 

Qscat is the scattering efficiency. Both quantities are computed using 

classical Mie theory. See G. F. Bohren, and D.R. Hoffman , Absorption and 

Scattering of Light by Small Particles , John Wiley, New York (1983) for 

20 additional information concerning this computation. 

The real and imaginary parts of 0(r s -r d ,go) can be used to predict the 

measured phase-shift, 0, and amplitude modulation, M, measured at detector 
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collected for the given X t and stored by processor 50. In one embodiment, 
phase error inherent in the instrumentation of system 20 is reduced by 
collecting the multiply scattered light in two locations, rdl and r<j2. and 
determining Orel =|0(i"dlHKrd2)l as more ful| y described in conjunction with 
5 stage 520. 

Next, conditional 516 is encountered. Conditional 516 tests whether 
frequency domain information for all wavelengths of interest have been 
collected. If not, the next wavelength ? vl is selected by appropriate 
adjustments to source 21 and control flows back to stage 512; closing loop 
io 510. After frequency domain information for all desired wavelengths has 
been collected, then conditional 516 is satisfied and control flows to stage 
520. 

In stage 520, the isotropic scattering coefficient for each of the 
wavelengths A., is determined. This determination is based on the diffusion 
1 5 equation model of light propagation. Under conditions of multiple scattering, 
the time-dependent light propagation can be accurately modeled by a 
diffusion equation written here in the frequency domain: 

-D(A)V 2 0(r.«) + 103 O(r.co)4- ja a (A.)O(r.0>) = 5(r - r s )m 

c n(^) (1) 

where 0(r, co) is the fluence of photons at position r within a propagating 

20 "photon density wave" modulated at frequency, co; Cn is the speed of light 

through the medium; m is the depth of sinusoidal modulation of the source 

located at position r s ; Ma(^) is the wavelength dependent absorption 

coefficient that is equal to£e[C], the sum of the product of extinction 
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j index to wavelengths X. 

M maximum number of wavelengths X. 

p discrete low pass filter index. 

5 VECTOR AND MATRIX VARIABLES 

I identity matrix, 

r position vector. 

r s position vector of the source fiber 30, [cm], 

rd position vector of the detector fiber 36 t [cm]. 

io C, Matrix containing updates to size distribution and volume 

fraction. 

£ Discrete filter parameter. 



Jacobian matrix describing changes in isotropic 
scattering with size distribution and volume fraction. 

15 

TABLE 1 Description of Selected Variables 

Referring to FIGS. 1 and 2, process 500 begins by entering loop 510 at 
stage 512. In stage 512, particles P in medium M are exposed to pulses of 
20 light of wavelength X i from source 21 via fiber 30; where "j" is the wavelength 
index form j=1, M. The location of the source light is designated by position 
vector r s . In stage 514, multiply scattered light resulting from the exposure in 
stage 512 is collected at detector location rd corresponding to the position of 
fiber 36. Frequency domain quantities relating to 9(rd) and M(rd) are 
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a regutarization parameter. 

minimization function. 
<|> dispersed phase or solids volume fraction. 

O photon fluence, equal to the local concentration of 

5 photons times the speed of light [photons/cm 2 sec]. 

0 phase-shift of amplitude modulated light relative to 

incident light, [degrees or radians]. 
6 re | relative phase-shift between two different detectors at 

different positions [degrees or radians]. 
10 X wavelength of light. 

[is scattering coefficient, [cm" 1 ]. 

ja a absorption coefficient, [cm' 1 ]. 

Hs'=(1-g)Hs isotropic scattering coefficient [cm -1 ]. 
a mean particle diameter. 

15 co modulation frequency. 

SUPERSCRIPT VARIABLES 
o observed, 
c computed. 



20 



SUBSCRIPT VARIABLES 
index to Jacobian bins, 
index to particle diameter. 
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• • 

GENERAL VARIABLES 
a,b ( c t cl parameters governing the size distribution and 

volume fraction for inverse solution. 

k wavelength number, k=2n//. [nm"1] 

5 Cn speed of light through random medium, [cm/sec]. 

D optical diffusion coefficient [cm]. 

f(x) particle size distribution (volume distribution) 

[dimensionless]. 
g mean cosine of angular scatter. 

10 r radius of particle, [microns]. 

x diameter of particle, [microns]. 

m depth of sinusoidal intensity modulated light. 

M amplitude modulation of re-emitted, multiply 

scattered light. 

15 n relative refractive index of particle. 

N number of Jacobian Matrix bins. 

ND number density of particles in suspension 

N* number of particle sizes for discrete filter application. 

q positive integer frequency multiple. 

20 q q is an intermediate variable, q = [(4TTn/A)(sin(0/2))]. 

Q scat scattering efficiency of single scatterer. 

q^ angular scattering efficiency. 

S(q) structure factor. 

z number of discrete particle size increments. 
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digitize the output of sensors 34, 38 for subsequent analysis by processor 50. 
Processor 50 may be an electronic circuit comprised of one or more 
components. Similarly, processor 50 may be comprised of digital circuitry, 
analog circuitry, or both. Also, processor 50 may be programmable, an 
5 integrated state machine, or a hybrid combination thereof. Preferably, 
processor 50 is comprised of one or more devices of the a digitally 
programmable variety. 

Processor 50 is coupled to at least one input device 54 and at least 
one output device 56. Preferably, input device 54 is a keyboard or input 
10 control of a conventional variety, and output device 56 is a Cathode Ray Tube 
(CRT) based video display, printer, or other image display system known to 
those skilled in the art. System 20 is well-suited to conveniently deliver and 
detect light pulses at a selected rate and wavelength in a laboratory setting. 
In other embodiments, the configuration of system 20 may vary to provide 
15 and sense light at selected wavelengths and with a desired time-varying 
intensity using techniques known to those skilled in the art. 

FIG. 2 is a flow chart depicting process 500 which is performed using 
system 20. Preferably, processor 50 is configured to execute a program that 
performs at least the more calculation intensive aspects of process 500. With 
20 regard to these calculations, selected variables employed in the description of 
process 500 are presented as follows in Table 1: 
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Where the source 21 pulse repetition rate is co and q is a whole number 
multiplier, heterodyne subsystem 40 gain modulates sensors 34, 38 at a 
harmonic qco plus a cross-correlation frequency offset, Aco, using conventional 
techniques. For one embodiment, an RF signal from frequency synthesizer 
5 44 (Marconi Instruments Signal Generator, model 2022A) is phase-locked to 
the pulse frequency of source 21 via divider 42 and outputs the gain 
modulation frequency of qco +Aco for amplification by poWer amp 46 (model 
1403LA, ENI f Rochester, NY). The output of power amp 46 modulates the 
signals from sensors 34, 38 resulting in a signal of frequency Aco from each 
io sensor 34, 38. Comparative frequency domain information, including phase- 
shift (6) and relative magnitude (M) (also called demodulation amplitude) for 
harmonic frequency qco of the pulsed laser light delivered to sensors 34, 38 is 
obtained by comparing these outputs. By changing q, measurement of 
phase shift (0) and relative magnitude (M) may be obtained as a function of 
15 modulation frequency. Furthermore, by varying the wavelength (A.) of the 
output from source 21 , the phase shift (0) and relative magnitude (M) may 
be obtained as a function of k. 

In one embodiment, the source 21 pulse repetition rate is 80 MHz 
which is divided by 8 by divider 42 and input to synthesizer 44 to provide the 
20 phase-locked loop. A 4 MHz pulse train is selected with picker 26 and the 
corresponding cross-correlation offset frequency is 100 Hz. Synthesizer 44 
is preferably a Marconi Instruments signal generator, model 2022A, and 
amplifier 46 is a model 1403LA available from ENI, Rochester, NY. 
Preferably, interface 152 is a conventional data acquisition module suitable to 
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A. and 0.5a. outputs are separated by a prism to provide two outputs 
collectively represented as beam 27. The X output is selected for 
wavelengths in the 720 to 900 nm range, and the 0.5A. output is selected for 
wavelengths in the 360 to 450 nm range. 



which splits beam 27 to direct about 80% along optical fiber 30 and the 
remaining 20% along fiber 32 to reference sensor 34. Dght directed along 
fiber 30 enters sample tank 60 to encounter particles P suspended in fluid 
medium M. Particles P are sufficiently concentrated in medium M to multiply 

io scatter light pulses entering by fiber 30. This multiple scattering is 

diagrammatically represented by scattering path SP. Multiply scattered light 
is collected by fiber 36, and is then directed to sample sensor 38. Preferably, 
sensors 34, 38 are Photomultiplier Tubes (PMTs), (Hamamatsu R928, 
Hamamatsu, Japan) and fibers 30, 32, 36 are 1000 mm optical fibers (HCP- 

15 M1000T-08, Spectran, Avon, CT); however, in other embodiments other 
sensor types and couplings are envisioned as would occur to one skilled in 
the art. 

Sensors 34, 38 are operatively coupled to heterodyne subsystem 40 
and processor 50 via interface 52. It should be appreciated that the light 
20 pulse train incident on sensors 34 and 38 may be described in terms of a 
number of intensity modulated frequencies which are multiples of the pulse 
repetition rate via Fourier analysis. Heterodyne subsystem 40 and processor 
50 provide frequency domain information about light collected from sensors 
34, 38 in terms of one or more of the intensity modulation frequencies. 



5 



Beam 27 from source 21 is directed to a glass slide beam splitter 28 
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rod and excite titanium ions of the rod. In response, a monochromatic, 
unidirectional optical pulse train of light pulses is produced which is 
designated beam 25. Preferably, beam 25 has generally equally spaced light 
pulses of 2 picoseconds (psec) Full Width Half Maximum (FWHM) at a 
5 repetition rate of about 80 Megahertz (MHz). The output power of beam 25 is 
preferably maintained at about 20% of the input power of beam 23. 

Beam 25 is directed to pulse picker 26. Picker 26 provides two 
functional capabilities: (1) to select an output pulse frequency less than or 
equal to the pulse repetition rate of pulses from input beam 25 and (2) to 
io selectively halve the wavelength (a.) of the input beam. The frequency 
selection capability is provided by directing beam 25 to an Acousto-Optic 
Modulator (AOM) crystal. The AOM crystal controllably diffracts pulses from 
the impinging beam 25 at an angle relative to the main beam in accordance 
with a electromagnetic Radio Frequency (RF) input. The rate of diffracted or 
1 5 "picked" pulses may be selected to provide a number of different output 
repetition rates less than or equal to the input frequency. 

Picker 26 also includes a Second Harmonic Generation (SHG) crystal 
to selectively produce output wavelengths at one half the input wavelength X. 
Picker 26 is configured so that beam 25 passes through the SHG crystal after 
20 the AOM crystal. The SHG crystal emits a single photon for every two 

photons absorbed. The energy level of the emitted photon is approximately 
double the energy of each of the two absorbed photons. As a result, the 
wavelength corresponding to the emitted photon is about one half and the 
frequency is about doubled compared to each of the two input photons. The 
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DESCRIPTION OF PREFERRED EMBODIMENTS 



For the purposes of promoting an understanding of the principles of the 
invention, reference will now be made to the embodiments illustrated in the 
drawings and specific language will be used to describe the same. It will 
5 nevertheless be understood that no limitation of the scope of the invention is 
thereby intended. Any alterations and further modifications in the described 
techniques, systems, and devices, and any further applications of the 
principles of the invention as described herein are contemplated as would 
normally occur to one skilled in the art to which the invention relates. 

10 FIG. 1 depicts an analysis system 20 of one embodiment of the 

present invention. System 20 includes a light source 21 which includes Argon 
ion laser 22, Titanium:Sapphire laser 24, and pulse conditioner 26. Laser 22 
generates a monochromatic output beam 23 using electrical discharges in a 
conventional manner. In one preferred embodiment, laser 22 is provided by a 

15 Beam-Lok model 2060 laser available from Spectra Physics. CA which is 
configured to deliver a beam with a diameter of less than 2 millimeters (mm) 
at a wavelength of 514.5 nanometers (nm). A constant output power of the 
Beam-Lok model 2060 laser is selected from a range of 0 and 10 watts (W) 
for this embodiment. 

20 Beam 23 from laser 22 is configured to "pump" laser 24 using 

conventional techniques. In one embodiment, source 21 includes a Tsunami 
Ti:Sapphire laser, model 3950B, available from Spectra Physics, Mountain 
View, CA for laser 24 and the Beam-Lok model 2060 for laser 22. Beam 23 
is directed into a cavity of laser 24 to impinge upon a titanium-doped sapphire 
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dispersion of Ti0 2 particles. 

FIGS. 12a-12b are graphs depicting size distributions for aqueous 
dispersions of Ti0 2 particles supplied as suspension samples S1 and S2, 
respectively, which compare measurement of size distribution in accordance 
5 with the present invention to measurements by Dynamic Light Scattering 
(DLS) and X-ray scattering techniques. 

FIG. 13 is a flow chart illustrating another process of the present 
invention capable of performance, for example, by the systems of FIGS. 1 
and 3. 

io FIG. 14 is a comparative schematic illustration of selected aspects of 

particle-to-particle interaction structure factor models of the present invention. 

FIG. 15 is a graph comparing selected structure factor models in terms 
of isotropic scattering coefficient and light wavelength. 

FIGS. 16a-16c are graphs comparing the relationship of the isotropic 
15 scattering coefficient with the percent solids by volume from polydisperse 
polystyrene particles in suspension having a mean diameter of about 0.4 
micron, 0.2 micron, and 0.15 micron, respectively. 

FIGS. 17a-17c are graphs depicting simulation results using a soft 
sphere particle-to-particle interaction structure factor model with simulated 
20 noise levels of 0%, 1%, and 5%, respectively. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 1 is a schematic view of an analysis system of one embodiment of 
the present invention. 

FIG. 2 is a flow chart illustrating one process of the present invention 
performed by the system of FIG. 1 . 

FIG. 2a is a graph illustrating a time domain measurement of a photon 
time-of-fiight distribution resulting from an incident light palse in a scattering 
media. 

FIG. 2b is a graph illustrating the frequency domain measurement of 
light propagation in a scattering media in terms of phase shift, 9, and 
amplitude, M. 

FIG. 3 is a schematic view of a chemical process control system of 
another embodiment of the present invention. 

FIG. 4 is a graph depicting the relationship between phase shift and 
modulation frequency for a dispersion of DOW PP722 particles. 

FIG. 5 is a graph depicting the isotropic scattering coefficient for a 
dispersion of DOW PP722 particles as a function of wavelength. 

FIGS. 6-8 are graphs depicting size distribution for DOW PP722, 
PP755, and PP788 particle dispersions without a mass balance constraint, 
respectively. 

FIGS. 9-1 1 are graphs depicting size distribution for DOW PP722, 
PP755, and PP788 particle dispersions with a mass balance constraint, 
respectively. 

FIGS. 12 is a graph depicting size distribution for an aqueous 
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algorithm to inversely solve for size distribution, volume fraction, or structure 
factor of a number of particles using observed values representative of the 
isotropic scattering coefficient of the particles at different wavelengths. This 
estimation may be based on an expected form of the distribution. Also, the 
5 estimation may account for mass of the particles. 

An additional object of the present invention is to provide a technique 
for determining a particle interaction parameter which accounts for particle-to- 
particle interactions that vary with particle concentration and influence light 
scattering by the particles. 
10 Further objects, advantages, features, forms, and aspects of the 

present invention will become apparent from the drawings and description 
contained herein. 
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provides a detected light signal corresponding to multiply scattered light from 
the particles at different wavelengths in response to exposure to the source. 
The processor is responsive to the emission signal and includes a calculation 
means for generating an output signal corresponding to at least one of a size 
5 distribution, volume fraction, or a structure factor of the particles in 

accordance with an observed isotropic scattering coefficient of the medium. 
The processor determines a value representative of the observed isotropic 
scattering coefficient from the detected light signal. The output device 
responds to the output signal to provide an output corresponding to at least 
io one of the size distribution, the volume fraction, or the structure factor of the 
particles. 

Accordingly, one object of the present invention is to determine particle 
size distribution, particle volume fraction, or a particle interaction parameter 
from multiply scattered light. 

15 Another object of the present invention is to provide a self-calibrating 

system for the on-line determination of particle size distribution and volume 
fraction by detecting the propagation characteristics of multiply scattered light 
from particles dispersed in a fluid medium. 

Yet another object is to determine one or more optical coefficients of a 

20 number of particles suspended in a medium by sensing propagation 

characteristics of multiply scattered light from the particles and to provide "a 
size distribution or volume fraction of the particles from one or more optical 
coefficients. 

Still another object of the present invention is to provide an estimation 
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expose the medium to a number of different predetermined wavelengths of 
incident light each having a predetermined time-varying intensity and a 
sensor spaced apart from the source. The sensor is configured to provide a 
detection signal corresponding to detected light. This detected light originates 
5 from the source and is multiply scattered by the particles. The system also 
includes a processor responsive to the first detection signal. The processor 
receives an exposure signal corresponding to the incident light and generates 
propagation signals by comparing the detection signal and the exposure 
signal for each of the wavelengths. These propagation signals characterize 
10 time of flight of the detected light through the medium resulting from multiple 
scattering by the particles. The processor is also configured to generate 
scattering signals corresponding to an isotropic scattering coefficient of the 
medium. The scattering signals are correspondingly determined from the 
propagation signals. The processor further generates an output signal 
i 5 indicative of one of size distribution or volume fraction of the particles which is 
determined from the scattering signals and a structure factor. This structure 
factor accounts for particle-to-particle interactions that influence light 
scattering behavior of the particles. The system further includes an output 
device responsive to the output signal to provide an output corresponding to 
20 the size distribution or volume fraction of the particles. 

Yet another form of the present invention is a system for analyzing a 
number of particles suspended in a medium in sufficient concentration to 
multiply scatter light. This system includes a light source, a sensor spaced 
apart from the light source, a processor, and an output device. The sensor 
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each other to multiply scatter light. Multiply scattered light is detected from 
the particles in response to the incident light to determine a first value 
corresponding to an observed isotropic scattering coefficient of the particles. 
An estimate corresponding to at least one of volume fraction or size 
5 distribution of the particles is established, and a second value is calculated as 
a function of the estimate. The second value corresponds to a calculated 
isotropic scattering coefficient. The first and second values are compared to 
establish an error, and the estimate is changed. The calculation of the 
second value, comparison of the first and second values, and changing of the 

10 estimate is repeated until the error reaches a desired minimum. 

In still another form of the present invention, a method of analysis 
includes exposing a fluid with a number of suspended particles to an incident 
light The particles are sufficiently concentrated in the fluid to scatter light. 
Multiply scattered light in response to this exposure is detected to determine a 

15 time-based value characteristic of propagation time of the multiply scattered 
light through the fluid. A quantity is determined as a function of the value 
which corresponds to an isotropic scattering coefficient of the fluid. An output 
is provided that corresponds to at least one of particle volume fraction, 
particle size distribution, or a particle interaction parameter corresponding to 
20 particle-to-particle interactions that influence light scattering behavior of the 
particles. The output is determined from the quantity. 

A further form of the present invention is a system for analyzing a 
number of particles suspended in a medium in sufficient concentration to 
multiply scatter light. This system includes a light source configured to 
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number of light wavelengths each being intensity-modulated at a 
predetermined frequency. Multiply scattered light from the medium in 
response to this exposure is detected to characterize propagation of the 
multiply scattered light through the medium with a number of values. The 
5 values each correspond to a different one of the wavelengths and are each 
representative of at least one of a phase or an amplitude of the multiply 
scattered light relative to the predetermined frequency. -An output is provided 
that is determined from the values. The output corresponds to at least one of 
particle size distribution, volume fraction, or a particle interaction parameter. 
10 This particle interaction parameter corresponds to a nonlinear relationship 
between particle concentration and an isotropic scattering coefficient for the 
particles. 

In an additional form, a number of particles suspended in a medium 
are exposed to a number of light wavelengths of predetermined time-varying 

15 intensity. Multiply scattered light form the particles in response to this 

exposure is detected to determine a number of values^each corresponding to 
a different one of the wavelengths. These values are each representative of 
a time of flight characteristic of the particles. The output corresponds to at 
least one of particle size distribution, volume fraction, or a particle interaction 

20 parameter. This particle interaction parameter corresponds to a nonlinear 
relationship between particle concentration and an isotropic scattering 
coefficient for the particles. 

In a further form, particles are exposed to an incident light with a 
predetermined time-varying intensity. The particles are sufficiently close to 
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More preferably, ^ is at least 100 times greater than ji a . Also, as used 
herein "multiply scattered light" refers to light that has traveled at least 5 times 
the mean isotropic scattering length for a particular wavelength, defined as 

5 When light is multiply scattered by particles in a fluid, it typically travels 

a greater distance and therefore has a longer travel time as compared to a 
direct path through the fluid without encountering the particles. The "time of 
flight" of photons through a scattering media is on the order of about 1 
nanosecond. Because many different scattering paths are likely, time of flight 
io is usually characterized as a distribution; however, this distribution typically 
varies with the light wavelength, the refractive character of the particles 
relative to the media, the size of the particles, and the volume fraction of the 
particles. Thus, in a further form of the present invention, measurements 
corresponding to photon migration time or travel time through the scattering 
15 media may be obtained and utilized to characterize the particles. 

At high particle concentrations, interactions between particles have a 
substantial influence on the light scattering properties of a multiple scattering 
media. This influence may degrade measurements which fail to account for 
particle-to-particle interactions. To adjust for particle interactions, a particle 
20 interaction parameter, which varies with particle concentration, may be 
utilized. In one form of the present invention, a particle structure factor is 
utilized, such as the Percus-Yevick model, which further provides an 
assessment of the stability of a corresponding particle suspension. 

Jn another form, particles suspended in a medium are exposed to a 
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SUMMARY OF THE INVENTION 
The invention relates to analysis of particles using photon migration 
techniques. Several aspects of the invention are novel, non-obvious, and 
provide various advantages. While the actual nature of the invention covered 
5 herein can only be determined with reference to the claims appended hereto, 
certain features which are characteristic of the present invention are 
described briefly as follows. 

In one form of the present invention, particle size distribution or volume 
fraction of a number of particles in suspension is determined from the 
10 measurement of one or more propagation characteristics of multiply scattered 
light re-emitted from the suspension. As used herein, "particles" include 
crystals, solids, liquid droplets, microbes, or microbe organelles. Also, as 
used herein, "multiple scattering" of light refers to the condition where, at a 
given wavelength, a photon is more likely to be consecutively scattered many 
is times between particles due to a refractive index mismatch between the 

dispersed particles and the host medium before being absorbed or detected. 
Generally, the multiple scattering condition entails that the particles be close 
enough to one another so that the distance a photon travels between 
scattering events is much smaller than the length traveled before absorption. 
20 Thus, multiple scattering of light occurs when n s '» |a a , where the absorption 
coefficient, n 3 , indicates the ability of a substance to absorb light of a 
particular wavelength, and the isotropic scattering coefficient, (1-g)m=^ s \ 
indicates the ability of a substance to scatter light of a particular wavelength in 
any particular direction. Preferably ^ is at least 10 times greater than n a . 
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various other system parameters. Also hampering the effectiveness of these 
techniques is the need to calibrate the equipment in situ. Process upsets, 
feedstock changes, or even normal batch process changes may invalidate the 
calibration. Depending upon the application, the sensor output for feedback 
5 control could be catastrophic absent proper calibration. 

Thus, a need remains for a technique to analyze particles in a process 
stream which is self-calibrating and interrogates the process stream without 
requiring sampling or dilution. Preferably, this technique may be used to 
determine particle size distribution and volume fraction regardless of particle 
10 concentration level. 
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The focal position is made to vibrate at a high rate so that the beam travels 
significantly faster than particles in the medium. In theory, particles intercept 
the beam and reflect light for a duration of time which is proportional to the 
particle size. The reflected light is detected and timed, and a particle size is 
5 inferred from this information. To reconstruct a size distribution, the 
instrument counts the number of times each size occurs. Besides the 
constraint on speed of the beam relative to particle motion, tests have shown 
that this system consistently measures small particles as too large and large 
particles as too small. Also, the range swept by the laser beam must be 
10 smaller than the smallest particle size to be measured. 

Similar to DLS and diffraction measurements, laser reflection 
techniques are based upon discrete single particle measurements to 
reconstruct PSD. Laser reflection measurements are sensitive to erroneous 
positioning of the focal plane, contribution of higher order scattering, and 
15 sampling error brought about by the hydrodynarnic partitioning of large 

particles away from the wall and sensor head. In addition, the measurements 
are reported not to be accurate with non-opaque particles or dispersed 
droplets. Sparks, R.G., and C.L. Dobbs, The Use of Laser Backscatter 
Instrumentation for On-Line Measurement of the Particle Size Distribution of 
20 Emulsions , Part. Part. Syst. Charact. 10, 279-289 (1993) provides additional 
information concerning this technique. 

Turbidity, DLS, Diffraction, and Laser reflection are all limited to some 
extent by the multiple scattering of light by the particles. These systems 
attempt to confine this problem by sampling and dilution, or by adjusting 
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statistical number of measurements can provide PSD's from diluted samples 
of particles. 

Turbidity, DLS, and diffraction measurements all require careful 
calibration of the light source and detector to provide meaningful 
5 measurements. Also, the possibility that wavelength dependent sample 
absorbance will vary during normal process disturbances and feedstock 
changes threatens the accuracy of these techniques. More importantly, these 
approaches all require a relatively dilute sample compared to the usual 
process concentration in order to be effective. 

10 Turbidity, DLS, and diffraction techniques suffer from other limitations 

which complicates on-line implementation. For example, on-line application 
of these techniques would require automated sampling procedures. For 
sizing solids, side stream measurements frequently create maintenance 
problems such as clogged pumps, conduits and filters within sampling 

15 devices. For sizing of liquid droplets, side stream measurements often induce 
coalescence and breakup. Furthermore the mechanical action of most 
automated sampling procedures may change the particle, crystal, or 
dispersed phase size distribution. Moreover, on-line measurements under 
these approaches may require substantial dilution, phase separation, or 
20 sample destruction to be effective. 

One system which attempts to solve these problems is a laser 
reflection or "backscatter" technique, such as the PAR TEC 100 available 
from Laser Sensor Technology. This system includes a laser with a narrow 
beam focused directly into a polydisperse medium undergoing processing. 
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and F.R. Hallet, Spherical Particle Size Determination by Analytical Inversion 
of the UV-Visible-NIR Extinction Spectrum , Appl. Optics, 35, 193-200 (1996); 
and J. Vavra, Antalik, J., and M. Liska, Application of Regression Analysis in 
Spectroturbidity Size Characterization Methods , Part. Part. Syst. Charact. 12, 
5 38-41 (1995) are cited as sources of additional information concerning 
turbidity based systems. 

Another optical approach is Dynamic Light Scattering (DLS), also 
termed quasi-elastic light scattering or photon correlation spectroscopy. DLS 
systems monitor the fluctuation of light intensity due to the Brownian motion 
10 of a single particle into and out of the near-field. From the time-dependence 
of intensity fluctuations, the particle diffusion coefficient can be computed and 
the radius obtained from the Stokes-Einstein equation. A statistical number of 
measurements of diluted and pretreated samples can provide a PSD. 
Notably, DLS techniques are often used as the laboratory "standard" for 
15 spherical particles having a diameter of less than 10 microns. Thomas, J.C. 
and V. Dimonie, Fiber Optic Dynamic Light Scattering from Concentrated 
Dispersions. 3: Particle Sizing in Concentrates , Appl. Optics., 29, 5332-5335 
(1990), and U.S. Patent Nos. 5,502,561 and 4,781,460 are cited as sources 
of additional background information concerning DLS techniques. 
20 In still another approach, angular light scattering or "diffraction" 

measurements are employed which monitor the angular scatter of light at a 
single wavelength due to diffraction from a single particle. Using classical 
scattering theory and known refractive indices of fluid and particle, an 
equivalent radius can be computed from an inverse solution. Again, a 
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ability to determine PSD and the volume occupied by the dispersed particles 
(the "volume fraction") through "on-line" observation of the particles as they 
participate in the process remains elusive. These limitations adversely 
impact the development of optimal process controls for batch, semi-batch, 
5 and continuous reactors. The absence of this on-line measurement capability 
also impedes accurate model verification important in the formulation of a 
control model for non-linear processes. Typically, the lack^pf accurate 
modeling coupled with the lack of robust measurements have forced the use 
of open-loop control based upon less efficient downstream measurements; 
10 adversely impacting quality and efficiency. 

Several laboratory particle sizing techniques exist, such as size- 
exclusion chromatography, capillary hydrodynamic fractionation, and 
photosedimentation; however, these techniques cannot provide on-line 
information about the process. Similarly, several optical techniques have 
15 been developed to provide PSD information. One of these techniques 
involves turbidity measurements which monitor the attenuation of light at 
multiple wavelengths traveling through a sample along a straight line path 
(180° relative to the incident light). Because these measurements do not 
account for backscatter nor multiple scattering of light back into the path 
20 length, turbidity analysis is only effective for diluted samples such that the 
product of turbidity, x, and pathlength, L, is less than about 0.3. In other 
words, the optical path of the sample is no less than about three times the 
mean distance between scattering particles. Van de Hulst, H. Light 
Scattering by Small Particles Dover Publications, New York (1983); J. Wang, 
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or another designated dimension. For example, Particle Size Distribution 
(PSD) is highly relevant to the application, texture, and appearance of 
titanium dioxide based paint products. In addition, emulsion polymerization 
processes produce paints, various coatings, and synthetic rubbers to name 
5 only a few. Since emulsion polymerization involves the growth of suspended 
polymer particles. PSD measurements can yield insight into the extent of 
reaction and molecular weight distribution and can also provide means for 
product characterization. These measurements can also be adapted to many 
crystallization processes, such as in the food industry, pharmaceutics, 
10 agricultural products, and bulk chemicals. One way to optimize many of the 
processes involving a dispersion of particles is by providing on-line 
measurement of particle size distribution in combination with robust process 
control responsive to these measurements. Rawlings, J.B., Miller, S.M., and 
W.R. Witkowski, Model Identification and Control of Solution Crystallization 
15 Processes , Ind. Eng. Chem. Res., 32, 1275-1296 (1993); Farrell, R.J. and Y- 
C. Tsai, Nonlinear Controller for Batch Crystallization: Development and 
Experimental Demonstration, AlChE J. t 41, 2318-2321 (1995); and 
Dimitratos, J. t Elicabe, G. t and C. Georgakis, Control of Emulsion 
Polymerization Reactors , AlChE J. t 40, 1993-2021 (1994) are cited as 
20 additional sources of background information concerning chemical process 
control. 

Despite the usefulness of PSD information, technology has not 
advanced such that these measurements can be made consistently under a 
wide variety of conditions and in a cost-effective manner. Moreover, the 
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CROSS-REFERENCES 
This application is a continuation-in-part of U.S Patent Application No. 
08/747,1 12, filed 8 November 1996 and claims the beneflt.of U.S. Provisional 
Application No. 60/050,809, filed 26 June 1997. 



The U.S. Government has a paid-up license in this invention and the 
right in limited circumstances to require the patent owner to license others on 
reasonable terms as provided by the terms of the National Science 
15 Foundation Young Investigator Award No. BES-9496239 and the National 
Institutes of Health Research Center Development Award K04-CA68374. 



20 scattered light, and more particularly, but not exclusively, relates to the 

determination of size distribution and volume fraction of particles using photon 
migration techniques. 

Today's chemical industry heavily relies on particulate or dispersed 
phase processes. The quality of many industrial products produced by these 

25 processes relates to the size distribution of the particles in terms of diameter 
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International application No. 



Date stamp of the receiving Office 



Applicant 



PURDUE RESEARCH FOUNDATION, et al. 



CALCULATION OF PRESCRIBED FEES 

I . TRANSMITTAL FEE 



240 



SEARCH FEE 
International search to be carried out by 



US 



450 



(If two or more International Searching Authorities are competent in relation to the international 
application, indicate the name of the Authority which is chosen to carry out the international search.) 



appli 

INTERNATIONAL FEE 

Basic Fee 102 
The international application contains 



. sheets- 



first 30 sheets 
72 



530 



10 



720 



remaining sheets additional amount 

Add amounts entered at b. and b, and enter total at B 



1250 



B 



Designation Fees 

The international application contains 
73 x 

number of designation fees 
payable (maximum 11) 



73 



128 



_ designations. 



1408 



amount of designation fee 



Add amounts entered at B and D and enter total at I 

(Applicants from certain States are entitled to a reduction of 75% of the 
international fee. Where the applicant is ( or all applicants are) so entitled, the 
total to be entered at I is 25% of the sum of the amounts entered at B and D.) 

k FEE FOR PRIORITY DOCUMENT 

5. TOTAL FEES PAYABLE 

,Add amounts entered at T. S, I and P. and enter total in the TOTAL box 



2658 



30 P 



$3,378.00 



TOTAL 



| | The designation fees are not paid at this time. 



MODE OF PAYMENT 

E authorization to charge 
deposit account (see below) 

GT1 cheque 

| [ postal money order 



| | bank draft 
| | cash 

| [ revenue stamps 



| | coupons 

[ | other (specify): 



DEPOSIT ACCOUNT AUTHORIZATION (this mode of payment may not be available at all receiving Offices) 
The RO/ US | [ is hereby authorized to charge the total fees indicated above to my deposit account. 

is hereby authorized to charge any deficiency or credit any overpayment in the total fees indicated above to my 
XI deposit account. 

is hereby authorized to charge the fee for preparation and transmittal of the priority document to the International 
Bureau of WIPO to my deposit account. 



□ 



23-3030- 



Deposit Account Number 



t account. s-\ 

ixs /n? *d S^So^^^ 



Date Idav/month/year) 



Signature L . Scott PAYNTER 



39 J _797_ 



Form I'CT/RO/MH (Annex ) (January reprint July 1 1 >97 ) 



Sec Xon>* to the fee calculation sheet 



w 


International Application No.^^ 


REOUEST 


International Filing Date 


The undersigned requests that the present 
international application oe processed 
according to the Patent Cooperation Treaty. 


Name of receiving Office and "PCT International Application" 




Applicant's or agent's file reference 

(if desired) (12 characters maximum) PUR-63/75P 


Box No. I TITLE OF INVENTION 

PARTICLE ANALYSIS SYSTEM AND METHOD 


Box No. n APPLICANT 


Name and address: (Famitv name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (i.e. country) of residence if no State of residence is indicated below.) 

PURDUE RESEARCH FOUNDATION 

Office of Technology Transfer 

1650 Engineering Administration Building 

Room 328 ENAD 

West Lafayette, Indiana 47906 
United States of America 


[ | This person is also inventor. 
Telephone No. 

Facsimile No. 

Teleprinter No ; 


State (i.e. country) of nationality: 

US 


State (i.e. country) of residence: 

US 


This person is applicant | 1 all designated rrrn 111 designated States except | I the United States 1 | the States indicated in 

for thr piirpn«« "f- I 1 Slates L^J ^ Unfed States of America | | of America only | | the Supplemental Box 


Box No. Ill FURTHER APPLICANT(S) AND/OR (FURTHER) INVENTOR(S) 


Name and address: (Family name followed by given name; for a legal entity, full official designation, 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant s State (i.e. country) of residence if no State of residence is indicated below.) 

SEVICK-MURACA, Eva 
7650 E . 100 N. 
Lafayette, Indiana 47905 
United States of America 


This person is: " 
| | applicant only 

|x | applicant and inventor 

1 1 inventor only (If this check-box 
' 1 is marked do not fill in below.) 


State (i.e. country) of nationality: 


State (i.e. country) of residence: 
US 


This person is applicant ( 1 all designated | | all designated States except tne United States | 1 the States indicated in 

for the purrns*"; r>f* 1 1 States 1 1 the United States of America UQfl of America only 1 1 the Supplemental Box 


pC | Further applicants and/or (further) inventors are indicated on a continuation sheet. 


Box No. IV AGENT OR COMMON REPRESENTATIVE; OR ADDRESS FOR CORRESPONDENCE 


The person identified below is hereby/has been appointed to act on behalf |~| agent | 1 common representative 
of the applicant(s) before the competent International Authorities as: UU ° \ — \ 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country.) 

PAYNTER, L . Scott 

W00DARD, EMHARDT, NAUGHTON, MORI ARTY & MCNETT 

Bank One Center/Tower, Suite 3700 

111 Monument Circle 

Indianapolis, Indiana 46204 US 

SEE CONTINUATION TO BOX NO. IV ON SHEET NO. 5 


Telephone No. 

317-634-3456 

Facsimile No. 

317-637-7561 

Teleprinter No. _ ^ 
810-341-3283 


1 — | Mark this check-box where no agent or common representative is/has been appointed and the space above is used instead to 
| 1 indicate a special address to which correspondence should be sent. . 



Form PCf/Re/101 (first sheet) (January 1997; reprint July 1997) See Notes to the request fc 



Continuation of Box No. Ill 



Sheet 



No. AS^^Ref 



If none 



k f HER APPLICANTS AND/OR (FUR 



f following sub-boxes is used* this sheet is not 



THEg) U 



INVENTORS 



eluded in the request. 



Name and address: (Family name * followed by given name: for a legal entity,^ 

The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (Le country) of residence if no State of residence is indicated below.) 

PIERCE, Joseph 
309 Juniper Street 
Lake Jackson, Texas 77566 
United States of America 



This person is: 

| | applicant only 

| x| applicant and inventor 

| | inventor only (If this check-box 
is marked, do not fill in below.) 



State (i.e. country) of nationality: 
US 



State (i.e. country) of residence: 
US 



This person is applicant 
for the purposes of: 



□ all designated I 1 all designated States except 
States I | the United States of America 



Hthe United States 
of America only 



□ the States indicated in 
the Supplemental Box 



Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (i.e. country) of residence if no State of residence is indicated below.) 

RICHTER, Steven 
16 Queens Court 
Brunswick, Georgia 31521 
United State's 65 America 


This person is: 

| | applicant only 

| X| applicant and inventor 

| | inventor only (If this check-box 
1 1 is marked do not fill in below.) 


State (i.e. country) of nationality: 
US 


State (i.e. country) of residence: 
US 



This person is applicant 
for the purposes of: 



□ all designated j | all designated States except 
States »r_ J o , * 



I | the United States of America 



r^-i the United States 
l-X I of America only 



□ the States indicated in 
the Supplemental Box 



Name and address: (Familv name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicants State (ie. country) of residence if no State of residence is indicated below.) 

SHINDE, Rajesh 
1901 Union Street #126 
Lafayette, Indiana * 47904 
United States of America 


This person is: 

| | applicant only 

|X | applicant and inventor 

1 1 inventor only (If this check-box 
■ 1 is marked, do not fill in below.) 


State (i.e. country) of nationality: 
IN 


State (Le. country) of residence; 
US 



This person is applicant | | all designated f— i all designated States except f— 1 the United States i 1 the States indicated in 

for the purposes of: I I States I 1 the United States of America I A l of America only | | the Supplemental Box 



Name and address: (Familv name followed by given name: for a legal entTty.full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (ie. country) of residence if ho State of residence is indicated below.) 

BALGI, Ganesh 
246 Longley Drive 
Lebanon, Indiana 46052 
United States of America 


This person is: 

| | applicant only 

|x | applicant and inventor 

| | inventor only (If this check-box 
b marked do not fill in below.) 


State (i.e. country) of nationality: 
IN 


State (Le. country) of residence: 
US 



This person is applicant 
for the purposes of: 



□ all designated r~ I all designated States except 
States | | the United States of America 



Sine United States 
of America only 



□ the States lndicaiea in 
the Supplemental Box 



y | Further applicants and/or (further) inventors are indicated on another continuation sheet. 



Form PCT/RO/ 101 (continuation sheet) (January 1997; reprint July 1997) 



See Notes to the request form 



Continuation of Box No. Ill 



Sheet No. 



APPLICANTS AND/OR (FURTHER) 



Agent's Ref : pur^63/75P 

c 



If none of the folloWH^ sub-boxes is used, this sheet is not to be inc< 



liumnn thei 



the request 



Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (i.e. country) of residence if no State of residence is indicated below.) 

* 

301 Huckleberry. Drive 
Lake Jackson, Texas 77566 
United States of America 


This person is: 

| | applicant only 

|X | applicant and inventor 

1 j inventor only (If this check-box 
1 1 is marked, do not fill in below.) 


State (Le. country) of nationality: 
US 


State (Le. country) of residence: 
US* 


This person is applicant 1 1 all designated 1 1 all designated States except | | th< 

for the purposes of: 1 1 States | | the United States of America LXl of 


1 United States | | the States indicated in 
America only | | the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity. Ml official destination. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (i.e. country) of residence if no State of residence is indicated below.) 

JIANG, Htiabei 

205-05 Airport Road 

West Lafayette, Indiana 47906 

United States of America 


This person is: 

| | applicant only 

|X | applicant and inventor 

j | inventor only (If this check-box 
is marked, do not fill in below.) 


State (i.e. country) of nationality: 

CN - 


State (i.e. country) of residence: 
US 


This person is applicant | | all designated ( 1 all designated States except f-71 the United States | | the States indicated in 

for the purposes of: l_l States | | the United States of America of America only | | the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (i.e. country) of residence if no State of residence is indicated below.) 


This person is: 

| [ applicant only 

| | applicant and inventor 

1 1 inventor only (If this check-box 
' ' is marked, do not fill in below.) 


State (i.e. country) of nationality: 


State (i.e. country) of residence: 


This person is applicant . 1 1 all designated I 1 all designated States except r — | the United States 1 1 the States indicated in 

for the purposes of: 1 1 States | | the United States of America LJ of America only | J the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity Jull official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicants State (Le. country) of residence if ho State of residence is indicated below.) 


This person is: 

| | applicant only 

| | applicant and inventor 

1 1 inventor only (If this check-box 
1 is marked do not fill in below.) 


Slate (ie. country) of nationality: 


State (Le. country) of residence: 


This person is applicant | 1 all designated | 1 all designated States except 1 — 1 the United States | 1 the S 

for the purposes of: 1 1 States | | the United States of America | | of America only | | the Supplemental Box 


[ | Further applicants and/or (further) inventors are indicated on another continuation sheet. 



Sheet No. 



Agent » a^Ref : PUR- 63/ 75P 



Box No.V 



DESIGNA 



The following designations 
Regional Patent 



TlflttOl 

is a^Wei 



F STATES 



lereby made under Rule 4.9(a) (mark the applicl 



ck- boxes; at least one must be marked)' 



0 

or 



AP 



EA 



ARIPO Patent: GH Ghana, KE Kenya, LS Lesotho, MW Malawi, SD Sudan. SZ Swaziland, UG Uganda, 
Z W Zimbabwe, and any other State which is a Contracting State of the Harare Protocol and of the PCT 
Eurasian Patent: AM Armenia, AZ Azerbaijan, BY Belarus, KG Kyrgyzstan. KZ Kazakstan. MD Republic of 
Moldova, RU Russian Federation. TJ Tajikistan, TM Turkmenistan, and any other State which is a Contracting State 
of the Eurasian Patent Convention and of the PCT 
IS EP European Patent: AT Austria, BE Belgium, CH and LI Switzerland and Liechtenstein, DE Germany, DK Denmark, 
ES Spain, FI Finland, FR France, GB United Kingdom, GR Greece, IE Ireland, IT Italy, LU Luxembourg, MC Monaco, 
NL Netherlands. PT Portugal, SE Sweden, and any other State which is a Contracting State of the European Patent 
Convention and of the PCT 

S OA OAPI Patent: BF Burkina Faso, B J Benin, CF Central African Republic, CG Congo, CI Cote d* Ivoire, CM Cameroon, 
GA Gabon. GN Guinea, ML Mali, MR Mauritania, NE Niger, SN Senegal. TD Chad. TG Togo, and any other State 
which is a member State of OAPI and a Contracting State of the PCT (if other kind of protection or treatment desired, specify 

on dotted line) 

National Patent (if other kind of protection or treatment desired, specify on dotted line): 



0 

a 
sr 

(3 

a 
ca 
ej 
a 

Ef 
BT 

® 
0 

H 

\s 

BT 
BT 

m 

BT 
IS 
SI 

a 
© 

a 
is 
a 
sr 

ar 
a 
a 
a 



AU 
AZ 
BA 
BB 
BG 
BR 
BY 
CA 



CU 
CZ 
DE 
DK 
EE 
ES 
FI 
GB 
GE 
GH 
HU 
IL 
IS 
JP 
KE 
KG 
KP 



AL Albania 

AM Armenia 

AT Austria 

Australia 

Azerbaijan 

Bosnia and Herzegovina 

Barbados 

Bulgaria 

Brazil 

Belarus 
Canada 

CH and LI Switzerland and Liechtenstein 

CN China 

Cuba 

Czech Republic 

Germany 

Denmark 

Estonia 

Spain 

Finland 

United Kingdom 

Georgia 

Ghana 

Hungary .* 

Israel 

Iceland 

Japan 

Kenya 

Kyrgyzstan 

Democratic People's Republic of Korea 



E 

a 

a 
ra 
m 
or 
0 
is 
m 

a 
eg 

El 

a 
sr 



LV Latvia 

MD Republic of Moldova 

MG Madagascar .„ 

MK The former Yugoslav Republic of Macedonia 

MN Mongolia 

MW Malawi 

MX Mexico 

NO Norway 

New Zealand 

Poland 

Portugal 

Romania 

Russian Federation 

Sudan 
Sweden 
Singapore 

Slovenia 

Slovakia 

Sierra Leone 

Tajikistan 

Turkmenistan 



NZ 
PL 
PT 
RO 
RU 
SD 
SE 
SG 
SI 
SK 
SL 
TJ 
TM 
TR 
TT 
UA 
UG 
US 



Turkey 

Trinidad and Tobago 

Ukraine /. 

Uganda 

United States of America . 



United States ot America 



KR 
KZ 
LC 
LK 
LR 
LS 
LT 
LU 



Republic of Korea 

Kazakstan 

Saint Lucia 
Sri Lanka 
Liberia 

Lesotho 

Lithuania 

Luxembourg 



EF UZ Uzbekistan 

VN Viet Nam 

gj YU Yugoslavia 

0T ZW Zimbabwe 

Check-boxes reserved for designating States (for the purposes of 
a national oatent) which have become party to the PCT after 
issuance of this sheet: 

Or ID Indonesia 



In addition to the designations made above, the applicant also makes under Rule 4.9(b) all designations which would be permitted 

under the PCT except the designation(s) of _ — : —— -r—\ 

The applicant declares that those additional designations are subject to confirmation and that any designation which is not contirmea 
before the expiration of 15 months from the priority date is to be regarded as withdrawn by the applicant at the expiration of that time 
limit. (Confirmation of a destination consists of the filing of a notice specifying that designation and the payment of the designation and confirmation 
fees. Confirmation must reuch the receiving Office within the / S-month time limit. ) 



Form PCT/RO/101 (second sheet) (July 1997) 



See Notes to the request form 



Sheet No. 



Agent 1 aHief : PUR-63/75P 



Supplemental Box If the Supple* 



ox is not used, this sheet need not be inciudt 



m 



e request. 



Use this box in the following cases: 

1. If in any of the Boxes, the space is insufficient to 
furnish all the information: 

in particular: 

(i) if more than two persons are involved as applicants 
and/or inventors and no "continuation sheet" is 
available: 



(ii) if. in Box No. II or in any of the sub-boxes of Box 
No. HI. the indication "the States indicated in the 
Supplemental Box" is checked: 



(Hi) if. in Box No. II or in any of the sub-boxes of Box 
No. III. the inventor or the inventor/applicant is not 
inventor for the purposes of alt designated States or 
for the purposes of the United States of America : 



( iv ) if in addition Ho the dgen t( s) indicated in Box No. I V, 
there are further agents: 

(v) if in Box No. V, the name of any State (or OA PI) is 
accompanied by the indication "patent of addition, " 
or "certificate of addition. " or if. in Box No. V, the 
name of the United States of America is accompanied 
by an indication "Continuation" or "Continuation- 
in-part": 

(vi) if there are more than three earlier applications 
whose priority is claimed: 



2. If the applicant claims, in respect of any designated 
Office, the benefits of provisions of the national law 
concerning non-prejudicial disclosures or exceptions to 
lack of novelty: 



in such case, write "Continuation of Box No — " [indicate the number 
of the Box] and furnish the information in the same manner as 
required according to the captions of the Box in which the space was 
insufficient; 

in such case, write "Continuation of Box No. Ill" and indicate for each 
additional person the same type of information as required in Box 
No. III. The country of the address indicated in this Box is the applicant 's 
State ( i.e. country) of residence if no State of residence is indicated below; 

in such case, write "Continuation of Box No. II" or "Continuation of Box 
No. Ill " or "Continuation of Boxes No. II and No. Ill" (as the case may 
be), indicate the name of the applicant s) involved and, next to (each) 
such name, the State(s) (and/on where applicable, ARIPO, Eurasian, 
European or OA PI patent) for the purposes of which the named person is 
applicant; 

in such case, write "Continuation of Box No. II" or "Continuation of Box 
No. Ill" or "Continuation of Boxes No. II and No. Ill" (as the case may 
be), indicate the name of the inventors) and, next to (each) such name, 
the State(s) (and/or, where applicable. ARIPO. Eurasian, European or 
OA PI patent) for the purposes of which the named person is inventor; 

in such case, write "Continuation of Box No. IV" and indicate for each 
further agent the same type of information as required in Box No. IV; 

in such case, write "Continuation of Box No. V" and the name of each 
State involved (or OA PI), and after the name of each such State (or 
OA PI), the number of the parent title or parent application and the date 
°f grant of the parent title or filing of the parent application; 



in such case, write "Continuation of Box No. VI" and indicate for each 
additional earlier application the same type of information as required 
in Box No. VI. 

in such case, write "Statement Concerning Non-Prejudicial Disclosures 
or Exceptions to Lack of Novelty" and furnish that statement below. 



Continuation to Box No. IV Agent 

WOODARD, Harold R.; EMHARDT, C. David; NAUGHTON, Joseph A., 
Jr.; MORI ARTY, John V.; McNETT, John C; HENRY , Thomas Q. ; 
DURLACHER, James M. ; REEVES, Charles R.; WAGNER, Vincent O.; 
ZLATOS, Steve; BEREVESKOS, Spiro; BAHRET, William F.; BROWNING, 
Clifford W;; FRISK, R. Randall; LUEDERS, Daniel J.; BECK, 
Michael D.; GANDY, Kenneth A.; THOMAS, Timothy N.; SISSELMAN, 
Kerry P.; JONES, Kurt N. ; ALL IE , John H. ; MICHAEL, Jeffrey A, ; 
BANTA, Holiday W. ; COLE, Troy J,; PAYNTER, L. Scott; LOWES, J. 
Andrew; HARRIS, Darrin Wesley; SCHANTZ, Matthew R . ; COY, 
Gregory B.; HIDAY, Lisa A.; DAN I LUCK, John V. and ROWE, James 
L., all of Woodard, Emhardt, Naughton, Moriarty & McNett, Bank 
One Center/Tower, Suite 3700, 111 Monument Circle, 
Indianapolis, Indiana 46204 United States of America 

l&ricJ Wo . ^/C5d, to 4 ] filed Z(o 3i*n<L. 



Form PCT7KO/10I (supplemental sheet) (January 1997; rcprim July 1997) 



See Notes to the request form 



Sheet No. 



Agent's Ref; PUR-63 /75P 

r 



Box No. VI PRIOR 



The priority of the followi 



IT^fl^VIM 

winHRier ap 



ier application(s) is hereby claimed: 



Further pri ority 1rr intlir atcd in the Supplemental Box [~j 
Tied* ~ 



Country 
(in which, or for which, the 
application was filed) 



Filing Date 
( day/month/year) 



Application No. 



Office of filing 
(only for regional or 
international application) 



item (I) 



US 



08 November 1996 
(08.11.96) 



08/747,112 



item (2) 



US 



26 June 1997 
(26.06.97) 



60/050,809 
(per postcard) 



item (3) 



Mark the following check-box if the certified copy of the earlier application is to be issued by the Office which for the purposes of the present international 
application is the receiving Office (a fee may be required): 

SThe receiving Office is hereby requested to prepare and transmit to the International , 1 . , _ x 
Bureau a certified copy of the earlier application(s) identified above as item(s) : (1J & (2) 



Box No. VII INTERNATIONAL SEARCHING AUTHORITY 



Choice of International Searching Authority (ISA) (If two or mine International Searching Authorities 

are competent to carry out the international search, indicate the Authority chttsen; the two-letter code may be used): ISA L 



US 



Earlier search Fill in where a search (international, international-type or other) by the International Searching Authority has already been carried 
out or requested and the Authority is now requested to base the international search, to the extent possible, on the results of that earlier search. Identify 
such search or request either by reference to the relevant application (or the translation thereof) or by reference to the search requesr. 
Country (or regional Office): ^ N g*|^ g§y&TM 



26 June 1997 (26.06.97) 



(per postcard) 



Box No. VIII CHECK LIST 



This international application contains 
the following number of sheets: 



1 . request 

2. description 

3. claims 

4. abstract 

5. drawings 

Total 



6 sheets 

57 sheets 

1 ^ sheets 

1 sheets 

22 sheets 



10 I. sheets 



This international application is accompanied by the item(s) marked below: 

5. |x | fee calculation sheet 



□ separate signed 
power of attorney 

□ copy of general 
power of attorney 

□ statement explaining 
lack of signature 

□ priority document(s) 
identified in Box No. VI 
as item(s): 



I 1 se P arat e indications concerning 

| | deposited microorganisms 

7 I I nucleotide and/or amino acid 
I I sequence listing (diskette) 

8- j5T] other (specify): Transmittal 
Letter (dup) 



Figure No. ,f, 



of the drawings (if any) should accompany the abstract when it is published. 



Box No. IX SIGNATURE OF APPLICANT OR AGENT 



Next to each signature, indicate the name of the person signing and the capacity in which the person signs (if such capacity is not obvious from reading the request). 

Applicant: Agent 
PURDUE RESEARCH FOUNDATION 
SEVICK-MURACA, Eva 

PIERCE, Joseph BALGI, Ganesh 

RICHTER, Steven KAO, Jeffery 

SHINDE, Rajesh JIANG, Huabei 





(L. Scott PAYNTER) 





Date of actual receipt of the purported 
international application: 


For rece 




2. Drawings: 


3. 


Corrected date of actual receipt due to later but 
timely received papers or drawings completing 
(he purported international application: 






| | received: 


4. 


Date of timely receipt of the required 
corrections under PCT Article 1 1(2): 






| | not received: 


5. 


International Searching Authority TC A . 
specified by the applicant: ior\ 1 




6. I 1 Transmittal of search copy delayed 
1— 1 until search fee is paid 





For International Bureau .use only 



Date of receipt of the record copy 
by the International Bureau: 



Form PCT/RO/I0I (last sheet) (January 1994: reprint July 1997) 



See Notes to the request form 



NOTIFICATION OF THE RECORDING 
OF A CHANGE 

(PCTRule92bis.1 and 
Administrative Instructions, Section 422) 



Pate of mai HngJday/mc 




Applicanfs or agent's file reference '^m^f^^ ^ 
PUR-63/75P v :;t . ¥ v.' •. 



International application No. 1 
PCT/US97/20539 



PAYNTER, L, Scott 
Woodard, Emhardt, Naughton, 
Moriarty & McNett .<* LV: 
Bank One Center/Tdwerv Suite 3700 ^ 
111 Monument Circle^ 

.Indianapolis, IN 4f>2l 

[EJATS-UNIS D'A^ERIQ^EI 



IMPORTANT NOTIFICATION , , 



international filing date (day/month/year) 

^ 07 November 1997L(0^jtl97)^ji 




Name and Address ^l^^^^^^^^ 

KAO, J effe ry ■,<, ^>^^ y ^^^-f^^:ifefe^:; s * ■ . ? 
301 Huckleberry Drive 
Lake Jackson, TX 77566 ■ ~^%<-v-; 



United States of America 




us 



State of Nationality T State of Residence L . 

'* f !mmmer» ™ jjwbw ******* 



Telephone No. - J - 








Facsimile No. 









Tele printer No. ^a***^.*^* ... „, .«^-^ 



2. The International Bureau hereby notifies the applicant that the following change hars been recorded concerning 
| | the person 



the name 



□ 



the address 



□ 



the nationality 



□ 



the residence ' 




Name and Address 

KAO, Jeffrey 
301 Huckleberry Drive ^ 
Lake Jackson, tX 77566-, 
United States of AmBring- 




State of Nationality 
US 



State of Residence 
US , . 



Telephone No. ^ . s:^-^^ 




Facsimile No. ^^^^^ :■- . ^ 



Teleprinter No. 



•••• -r<f^m^^.^>.- 



3. Further observations, if necessary:! 



4. A copy of this notification has : beehliVitf to^ 

Sthe receiving Office 
| J the International Searching Authority 

□ the International Prelimingry ^a^^hof^^^^ 



the designated Offices concerned 



□ 

| | the elected Offices concerned 
| | other 




Form PCT/IB/306 (March 1994) X:. . ? : f^^« 



v ^002006586^/ r 
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PA^JT COOPERATION TREAT^ 



From the INTERNATIONAL BUREAU 



PCT 

NOTIFICATION OF ELECTION 

(PCT Rule 61.2) 


To: 

United States Patentand Trademark 

Office 

(Box PCT) 

Crystal Plaza 2 

Washington, DC 20231 

ETATS-UNIS D'AMERIQUE 

in its capacity as elected Office 


Date of mailing (day/month/year) 
03 July 1998 (03.07.98) 




International application No. 
PCT/US97/20539 


Applicant's or agent's file reference 
PUR-63/75P 


International filing date (day/month/year) 

07 November 1997 (07.11.97) 


Priority date (day/month/year) 

08 November 1996 (08.1 1 .96) 


Applicant 

SEVICK-MURACA, Eva et al 



1. The designated Office is hereby notified of its election made: 

| X | in the demand filed with the International Preliminary Examining Authority on: 

03 June 1998 (03.06.98) 

| | in a notice effecting later election filed with the International Bureau on: 



2. The election 



□ 



was not 



made before the expiration of 19 months from the priority date or, where Rule 32 applies, within the time limit under 
Rule 32.2(b). 



The International Bureau of WIPO 


Authorized officer 


34, chemin des Colombettes 


S. De Michiel 


1211 Geneva 20, Switzerland 




Facsimile No.: (41-22) 740.14.35 


Telephone No.: (41-22) 338.83.38 



Form PCT/IB/331 (July 1992) 
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RECEIVED 

.JUL 1 3 1998 



Woodard. £mn 3 rd:. Naugmnn 



rT^NT COOPERATION TREM"^ 
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From the INTERNATIONAL BUREAU 



PCT 



INFORMATION CONCERNING ELECTED 
OFFICES NOTIFIED OF THEIR ELECTION 

(PCT Rule 61.3) 



Date of mailing (day/month/year) 
03 July 1998 (03.07.98) 



Applicant's or agent's file reference 
PUR-63/75P 



International application No. 
PCT/US97/20539 



To: 



PAYNTER, L., Scott 

Woodard, Emhardt, Naughton, 

Moriarty & McNett 

Bank One Center/Tower, Suite 3700 

111 Monument Circle 

Indianapolis, IN 46204 

ETATS-UNIS D'AMERIQUE 



IMPORTANT INFORMATION 



International filing date (day/month/year) 
07 November 1997 (07.11.97) 



Priority date (day/month/year) 

08 November 1996 (08.11.96) 



Applicant 



PURDUE RESEARCH FOUNDATION et al 



1. The applicant is hereby informed that the International Bureau has, according to Article 31(7), notified each of the following 
Offices of its election: 

AP :GH,KE,LS,MW,SD,SZ,UG,ZW 

EP :AT,BE,CH,DE,DK,ES,FI,FR,GB # GR-,rE,IT>LU,MC,NL,PT,SE 

National lAU^CBR^CA^CN.CZ^E^BJLJ^KP^KR^N^CNZ^PL^RCRU^^SK.US, 
VN 

2. The following Offices have waived the requirement for the notification of their election; the notification will be sent to them 
by the International Bureau only upon their request: 

EA :AM,AZ,BY,KG / KZ,MD,RU,TJ,TM 

OA :BF3J,CF,CG,CI,CM,GA,GN,ML,MR,NE,SN/TD,TG 

National ^L^M^AZ^BA.BB.BY^H.CU^K^E^ES^FLGE.GH^UJDJS^E^CKZ, 

LCXK,LR,LS,LTXUXV # MD,MG,MK,MW,MX,PT,SD,SG,SI,SL,TJ,TM,TR,TT,UA,UG,UZ, 

YU,ZW 

3. The applicant is reminded that he must enter the "national phase" before the expiration of 30 months from the priority date 
before each of the Offices listed above. This must be done by paying the national fee(s) and furnishing , if prescribed, a 
translation of the international application (Article 39(1 )(a)), as well as, where applicable, by furnishing a translation of any 
annexes of the international preliminary examination report (Article 36(3)(b) and Rule 74.1). 

Some offices have fixed time limits expiring later than the above-mentioned time limit. For detailed information about the 
applicable time limits and the acts to be performed upon entry into the national phase before a particular Office, see Volume I 
of the PCT Applicant's Guide. 

The entry into the European regional phase is postponed until 31 months from the priority date for all States designated for 
the purposes of obtaining a European patent. 



The International Bureau of WIPO 
34, chemin des Colombettes 
■ ~~ 1211 Geneva 20, Switzerland 


Authorized officer: 

S. Deyl 


vjichi&l 




Facsimile No. (41-22) 740.14.35 


Telephone No. (41-22) 


6.83.38 \ 




Form PCT/IB/332 (September 1997) 
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From the 

INTERNATIONAL PRELIMINARY 



PATENT COOPERATION TREA 

fcNING AUTHORITY 



Ta 

L. SCOTT PAYNTER 

WOODARD, EMHAROT, NAUGHTON. MORI ARTY 
BANK ONE CENTER/TOWER, SUITE 3700 
111 MONUMENT CIRCLE 
INDIANAPOLIS IN 46204 


TCT 

NOTIFICATION OF RECEIPT 
OF DEMAND 

(PCT Rule 61.1(b), first sentence 
and Administrative Instructions, Section 601) 


Date of mailing 

(daylmonthfyear) £2 vUN 1998 


Applicant's or agent's file reference 
PUR-63/75P 


IMPORTANT NOTIFICATION 


International application No. 

PCT/US97/20539 


International filing date (daylmonthfyear) 
07 NOV 97 


Priority date ( daylmonthfyear) 

08 NOV 96 


Applicant 

PURDUE RESEARCH FOUNDATION 



The applicant is hereby notified that this International Preliminary Examining Authority considers the following date as the 
date of receipt of the demand for international preliminary examination of the international application: 

This date of receipt is: 

the actual date of receipt of the demand. 
j j the date on which the proper corrections to the demand were timely received. 



3. | j This date is AFTER the expiration of 19 months from the priority date. 

Attention: The election(s) made in the demand does (do) not have the effect of postponing the commencement of the 
national phase until 30 months from the priority date (or later in some Offices) (Article 39( 1)). Therefore, the acts for 
entry into the national phase must be performed within 20 months from the priority date (or later in some Offices) 
(Article 22). 

For details, see Annex B to Form PCT/DB/301 sent by the International Bureau and Volume n of the PCT Applicant's 
Guide. 



r~| This notification confirms the information given in person or by telephone on: 



4. Only where paragraph 3 applies, a copy of this notification has been sent to the International Bureau. 



Name and mailing address of the EPEA/US 
Assistant Commissioner for Patents 
Box PCT 

Washington, D.C 20231 Attn: IPEA/US 

Facsimile No. 



Authorized officer 



Telephone No.(^7 Q3 ) 3Q S ~ 3b $ I 



Form PCT/IPEA/402 (July 1992) 
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PATENT COOPERATION TR EATY — 

rec • d 19 FEB 1999 



wipo 



PCT 

International preliminary examination report 

(PCT Article 36 and Rule 70) 



PCT 



Applicant's or agent's file reference 
PUR-63/75P 


FOR FURTHER ACTION See Notification of Transmittal of^t^tionsLX 

Preliminary Examination Report (Form PCT/IPEA/4 16) 


International application No. 
PCT/US97/20539 


International filing date (day/month/year) 
07 NOVEMBER 1997 


Priority date (day /month/year) 
08 NOVEMBER 1996 


International Patent Classification (IPC) or national classification and IPC 
IPC (6): GO IN 21/51 and US CI.: 356/336 


Applicant 

PURDUE RESEARCH FOUNDATION 



1. This international preliminary examination report has been prepared by this International Preliminary 
Examining Authority and is transmitted tcythe applicant according to Article 36. 



2. This REI 



' consists of a total of 



sheets. 



'This report is also accompanied by ANNEXES, i.e., sheets of the description, claims and/or drawings which have 
been amended and are the basis for this report and/or sheets containing rectifications made before this Authority, 
(see Rule 70. 16 and Section 607 of the AcUninistrative Instructions under the PCT). 

These annexes consist of a total of ^ S heets. 



3. This report contains indications relating to the following items: 
Basis of the report 
Priority 



I 




II 


□ 


III 


□ 


IV 


□ 


V 




VI 


□ 


VII 


□ 


VIII 


□ 



citations and explanations supporting such statement 



Date of submission of the demand 
03 JUNE 1998 


Date of completion of this report 
21 JANUARY 1999 


Name and mailing address of the IPEA/US 

Commissioner of Patents and Trademarks 
Box PCT 

Washington. D.C. 20231 
Facsimile No. (703) 305-3230 


Authorized officer , ^^-^ "> 

RICHARD ROSENBE^ofc^/'T — A 
Telephone No. (703) 308^u£*_ * J S \J 



Form PCT/IPE A/409 (cover sheet) (January 1994) 
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INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



International application No. 
PCT/US97/20539 



I. Basis of the regort 



1 . This report has been drawn on the basis of (Substitute sheets vrfuch have been furnished to the receiving Office in response to an invitation 
under Article 14 are referred to in this report as "originally filed" and are not annexed to the report since they do not contain amendments): 



| | the international application as originally filed. 

[~x| the description, pages (See Attached) , as originally filed. 



0 



pages 
pages 
pages 



, filed with the demand. 
, filed with the letter of . 
, filed with the letter of _ 



the claims, 



the drawings, 



Nos. (See Attached) 

Nos. 

Nos. 

Nos. 

Nos. 

sheets/6§ (See Attached) 

sheets/ftg 

sheets/ftg 

sheets/ftg 



, as originally filed. 
, as amended under Article 19. 
, filed with the demand. 

, filed with the letter of 

, filed with the letter of 



, as originally filed. 
, filed with the demand. 
, filed with the letter of, 
, filed with the letter of , 



2. The amendments have resulted in the cancellation of: 

|~x] the description, pages none 

|~x] the claims, 



Nos. 



36 



[~x] the drawings, sheets/ftg "one . 

3 . | I This report has been established as if (some of) the amendments had not been made, since they have been considered 
— to go beyond the disclosure as filed, as indicated in the Suppl e m e ntal Box Additional observations below (Rule 70.2(c)). 



4. Additional observations, if necessary: 
NONE 
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V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 



1. STATEMENT 
Novelty (N) 

Inventive Step (IS) 
Industrial Applicability (IA) 



Claims 1-35, 37-50 



Claims none 



Claims 1-35, 37-50 



Claims none 



Claims 1-35, 37-50 



Claims none 



YES 
NO 

YES 
NO 



YES 
NO 



2. CITATIONS AND EXPLANATIONS 

Claims 1-35, 37-50 meet the criteria set out in PCT Article 33(2)-(4), because the prior art does not teach or fairly 
suggest what is claimed. Although Paithankar et a! teaches in general the use of multiple wavelength scattering measurement 
for determining size distributions, there is no disclosure therein relating to specific details of how to obtain such 
measurements from such measurements and of specific claimed details for making such measurements. For example, the 
reference does not teach a method in which isotropic scattering coefficient is calculated as a function of an estimate of the 
volume fraction and/or the size distribution (claims 1-7, 32-35, 37-40 and 46-50). The reference does not discuss the 
generation of the "comparisons signals" (claims 8-13). The reference does not disclose the making the calculations 
iteratively (claims 1-7, 14-18, and 22-24), nor the use of a P-Y (Perkus-Yervick) hard sphere model (claims 17, 25, 30-31, 
40, 44, 48) nor of a Weibull distribution (claims 7, 18). The system of the reference does not disclose the use of a number 
of wavelengths "each being intensity -modulated at a predetermined frequency" and the detection signals "each being 
representative of at least one of a phase or an amplitude . . . relative to the predetermined frequency" (claims 19-25, and 
37-39). The reference does not teach "determining the particle interaction parameter" (claims 13, 26-31 and 34-35), or the 
calculation of the isotropic scattering coefficient as a function of an estimate of the volume fraction and the size distribution 
(claim 36). The reference does not teach the use of a structure factor in the calculation (claims 41-45). 



NEW CITATIONS 



NONE 



Form PCT/IPEA/409 (Box V) (January 1994)* 




This Page Blank (uspto) 



INTERNATIONAL PRELIMINARY EXAMINATION REPORT 
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Supplemental Box 

(To be u sed when tho-space in any of the preceding boxes is not sufficient) 

Continuation of: Boxes I - VIII Sheet 10 

I. BASIS OF REPORT: 

This report has been drawn on the basis of the description, 
pages, 1-29, 31-33, 36-38, 40-57, as originally filed, 
pages, none, filed with the demand, 
and additional amendments: 

pages 30, 34, 35, 39, filed with the letter of 22 January 1999. 

This report has been drawn on the basis of the claims, 

numbers, none, as originally filed. 

numbers, none, as amended under Article 19. 

numbers, none , filed with the demand. 

and additional amendments: 

Claims 1-35, 37-50, filed with the letter of 22 January 1999. 

This report has been drawn on the basis of the drawings, 

sheets, 1-22, as originally filed. 

sheets, none, filed with the demand. 

and additional amendments: 

NONE 
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where it is desired to minimize x 2 so that it is less than or equal to an 
established error amount. Conditional 546 tests whether this minimization 
has taken place. If the minimum has been reached, then control flows to 
stage 550 and the estimate is provided as the measured product <|>f(x). 
5 Subsequently, f(x) and may be separated by performing a summation 

CO 

approximation of |(()f(x)dx = <|> . 

o 

However, if conditional 546 is not satisfied, then control flows to stage 
548 to improve the estimate of §f(x) . One embodiment for rapidly converging 
to the desired minimum error is through application of a Jacobian matrix 
io relationship. Specifically, the product .4>f(x) is divided into N number of bins 
representing <j>f(x) h . where h = 1 to N. In order to update <|>f(x) h in each bin, 
the Jacobian matrix which describes the sensitivity of the isotropic scattering 
coefficient to changes in <t>f(x) h in each bin h and at each wavelength A., is 
computed. This Jacobian matrix is given by: 

| S(<|)f(x)) 1 <5((K(x)), A <5(<j>f(x)) N 

e ^X 

3 = j 3(4*00), a(4>f(x)) 2 A c(<f>f0O) N 
1 d(»% ° 6^% 

and each element is computed through a summation approximation of ~~ 
equation (3) with <f>f(x) h and <|>f(x) h + A<)>f(x) h . The updates, 5<|)f(x) h , 
resulting from the differences between measured and computed isotropic 
scattering coefficients, can then be obtained from: 
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10 



assuming an expected form of the PSD. Generally, particulate processes can 
usually be characterized by Gaussian or log-normal distributions, both of 
which can be represented by proper parametric choices for a Weibull 
distribution described as follows: 



c r (x-a) T- r rx-ayi 



(15) 



where a, b and c describe the peak location, width, and shape, respectively, 
of the distribution. For this embodiment, d corresponds to <j> , and the 
equation (3) summation approximation becomes: 



d- 



exp - 



(16). 



where the expected range of particle diameter x from x min to x max is determined 
and incremented into k number of discrete sizes Xj = x min + (i-t) Ax, 
Ax ='(x„,ax - x min ) / (z - 1) for i=1 to z, and the terms Q sca A jt n ' x i) and g^.n.x,) 
15 are determined from look-up tables based on k v n, x } . 

Using Newton's method, the updates based estimation of the a,b,c,d 
parameters can be obtained from the resulting system of equations: 

3 r 3A<; = o'k^) 0 -^;) 6 ] ^ 17 j _ 

where the Jacobian matrix 3 now represents the sensitivity of isotropic 
20 scattering coefficients measured at wavelengths j=1,M upon the four 
parameters (a, b, c and d): 
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facKK a(n;)j a(np< aai;)< 1 



1 da 


db 


dc 


ad 


1 








1 5a 


db 


dc 


ad 


WX 








L da 


db 


5c 


ad 



^ \ ua uu c/c uu 

I 



(18) 

The elements of the Jacobian are numerically computed using 
equation (3). The vectors, (Ms )° and (Ms ) c contain the observed and 
computed values of isotropic scattering coefficients at the wavelengths of 
5 interest: 

(K) n = ((^)^(rt) € ^-»(^)i M ) ( (K) c -((K)v(K)L.-(^)ij 

(19) 

AC, is the vector updating the four parameters, a f b, c and d: 
A£ = (5a,6b,5c,5d) T (20) 
io Thus, the particle sizing task now becomes to recover the four-parameters (a, 
b, c and d) to describe f(x) and <j> using the approach described in connection 
with loop 540. 

In still another embodiment, a low pass filter is employed to smooth the 
estimate based on the parameters for an expected form of distribution, such 
15 as the Weibull distribution. In this embodiment, the size distribution is 

computed from the parameters and then subjected to a digital low pass filter 
by averaging over a window of N*Ax. 

f(x p r w = a -^)f(x p r ,d z 2 fa,) 01 " 

,=P ~T (21) 
where ^ is a factor between 0 and 1 , and the summation is over the N* 
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intensity. In other embodiments of the present invention, the light source 
arrangement is configured to provide wavelengths and time-varying intensity 
suitable for interrogation of the particular process being monitored using 
techniques known to those skilled in the art. 
5 Analyzer 150 also includes optical fiber 1 34 to detect re-emitted light 

from the particles in vessel 160 after multiple scattering. It is preferred that 
analyzer 150 be of rugged design suitable for industrial applications. 
Analyzer 150 includes a processor configured to execute software to 
determine values representative of PSD or volume fraction in accordance 

io with the present invention. Furthermore, analyzer 150 is also operatively 
coupled to control elements 170 which are used to regulate processing in 
vessel 160. By way of non-limiting example, elements 170 may include 
valves, heaters, or agitation devices electronically regulated by analyzer 150. 
Analyzer 150 is further configured with appropriate programming and 

1 5 interfaces to provide one or more output signals to control elements 170 as a 
function of PSD or volume fraction. Thus, system 120 provides a closed loop 
feedback control system capability for regulating chemical processes from on- 
line measurement of PSD and/or volume fraction. 

Systems 20 and 120 of FIGS. 1 and 3, respectively, are described with 

20 corresponding light source fibers 30, 130 and detecting fibers 36, 136 

geometrically arranged to approximate infinite boundary conditions -- such 
that light travelling from fibers 30, 130 to fibers 36, 136 does not encounter a 
substantial boundary. Equations (6) or (7) are applied to determine the 
isotropic scattering coefficient p 5 ' and the absorption p a under such infinite 
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CLAIMS 

What is claimed is: 

1. A method of particle analysis, comprising: 

(a) exposing a number of particles in suspension to an incident light 
5 with a predetermined time-varying intensity, the particles being sufficiently 

close to one another to multiply scatter light; 

(b) detecting multiply scattered light from the particles in response to 
the incident light to determine a first value corresponding to an observed 
isotropic scattering coefficient of the particles; 

10 (c) establishing an estimate corresponding to at least one of volume 

fraction or size distribution of the particles; 

(d) calculating a second value from the estimate, the second value 
corresponding to a calculated isotropic scattering coefficient; 

(e) comparing the first and second values to establish an error; 
is (f) changing the estimate; and 

(g) repeating said calculating, comparing, and changing until the error 
reaches a desired minimum. 



2. The method of claim 1 , wherein the incident light includes a 

20 number of different wavelengths of light, the first value is determined for each 
of the wavelengths, and said establishing, said calculating, said comparing, 
and said changing are performed for each of the different wavelengths. 

3. The method of claim 1, further comprising controlling a process 
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as a function of the estimate and determining the estimate to generally 
maintain mass balance. 

4. The method of claim 1 , wherein the estimate corresponds to 
5 one of a number of parameters for an expected form of the size distribution 



5. The method of claim 4, wherein the incident light includes a 
number of different wavelengths of light, the first value is determined for each 
of the wavelengths, the parameters correspond to a Weibull distribution, and 
io said calculating includes: 

(i) establishing a number of particle size increments; 

(ii) determining a scattering efficiency and a mean cosine of 
scattering angle as a function of the wavelengths and the increments; and 

(iii) performing a first summation over a range of the particle sizes 
15 for each of the wavelengths, the first summation having the parameters, the 

scattering efficiency, and the mean cosine of scattering angle as arguments. 



6. The method of claim 5, wherein said comparing includes 
performing a second summation over a range of the wavelengths, the second 

20 summation including a numerical difference between the first and second 
values as an argument. 

7. The method of claim 5, wherein said changing includes updating 
the parameters as a function of the wavelengths, the first and second values, 
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and a third value corresponding to mass of the particles. 

8. A system for analyzing a number of particles suspended in a 
medium in sufficient concentration to multiply scatter light comprising: 
5 a light source with a predetermined time-varying intensity configured to 

expose the medium to a number of different predetermined wavelengths of 
incident light; 

a first sensor spaced apart from said source, said first sensor being 
configured to provide a first detection signal corresponding to detected light, 

io the detected light being multiply scattered by the particles; 

a processor responsive to the first detection signal and being 
configured to receive an exposure signal corresponding to said incident light, 
said processor being configured to generate: (a) a number of comparison 
signals each corresponding to a difference with respect to time between the 

15 detected light and the incident light for a corresponding one of the 

wavelengths, (b) a number of scattering signals each correspondingly 
determined from the comparison signals and each corresponding to an 
observed isotropic scattering coefficient of the medium for a different one of 
the wavelengths, and (c) an output signal indicative of one of size distribution 

20 or volume fraction of the particles, said output signal being determined as a 
function of said scattering signals; and 

an output device responsive to said output signal to provide an output 
corresponding to the size distribution or volume fraction of the particles. 
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9. The system of claim 8, further comprising a reaction vessel 
containing the particles and medium, and wherein said output device 
includes a control element operatively coupled to said reaction vessel and 
responsive to said output to regulate a process involving the particles. 

5 

10. The system of claim 8, wherein said processor is further 
configured to determine said size distribution or volume fraction as a function 
of mass of the particles. 

io 11. The system of claim 8, further comprising a second sensor 

providing a second detection signal, said comparison signals being 
determined as a function of said first and second detection signals for each of 
the wavelengths. 

15 12. The system of claim 8, wherein said processor further 

determines a structure factor indicative of particle-to-particle interactions, 
said structure factor varying in accordance with concentration of the particles 
in the medium. 

20 13. The system of claim 8, wherein said processor is configured to 

generate a particle interaction signal representative of particle-to-particle 
interactions that vary with particle concentration and influence light scattering 
behavior of the particles. 
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14. A system for analyzing a number of particles suspended in a 
medium in sufficient concentration to multiply scatter light, comprising: 

(a) a light source with a predetermined time-varying intensity 
configured to expose the medium to a number of different wavelengths of 

5 light; 

(b) a sensor spaced apart from said source, said sensor being 
configured to provide a detected light signal corresponding to multiply 
scattered light from the particles at the different wavelengths in response to 
exposure to said source; 

io (c) a processor responsive to said detected light signal, said processor 

including a calculation means for generating an output signal corresponding 
to at least one of particle size distribution, particle volume fraction, or a 
structure factor determined in accordance with an observed isotropic 
scattering coefficient of the medium, said structure factor being representative 

15 of particle interactions that influence light scattering behavior of said particles, 
said processor being configured to determine a value representative of said 
observed isotropic scattering coefficient from said detected light signal; and 

(d) an output device responsive to said output signal to provide an 
output corresponding to said at least one of said particle size distribution, said 

20 particle volume fraction, or said structure factor of the particles. 

15. The system of claim 14, further comprising a reaction vessel 
containing the particles and the medium, and wherein said output device 
includes a control element operatively coupled to said reaction vessel and 
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responsive to said output to regulate a reaction involving the particles. 

16. The system of claim 14, wherein said calculation means 
includes an estimating means for iteratively determining said size distribution 

5 or said volume fraction as a function of mass of the particles. 

17. The system of claim 14, wherein said structure factor is 
dependent on said particle size distribution and said particle volume fraction, 
and corresponds to a P-Y hard sphere model. 

18. The system of claim 14, wherein said calculation means 
includes a number of parameters corresponding to a Weibull distribution. 

19. A method of particle analysis, comprising: 

(a) exposing a medium with a number of suspended particles to a 
number of light wavelengths, the wavelengths each being intensity-modulated 
at a predetermined frequency; 

(b) detecting multiply scattered light from the medium in response to 
said exposing to characterize propagation of the multiply scattered light 
through the medium with a number of values, the values each corresponding 
to a different one of the wavelengths and each being representative of at 
least one of a phase or an amplitude of the multiply scattered light relative to 
the predetermined frequency; and 

(c) providing an output determined from the values, the output 
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corresponding to at least one of a particle size distribution, particle volume 
fraction, or a particle interaction parameter, the particle interaction parameter 
corresponding to a nonlinear relationship between particle concentration and 
an isotropic scattering coefficient for the particles. 



20. The method of claim 19, further comprising controlling a process 
in accordance with the output. 

21. The method of claim 19, wherein said providing includes 

io determining the observed isotropic scattering coefficient and an absorption 
coefficient for the particles for each of the wavelengths from a corresponding 
one of the values. 



particles for each of the wavelengths from a corresponding one of the values; 

(ii) establishing an estimate corresponding to at least one of the 
size distribution or the volume fraction; 

(iii) determining a calculated isotropic scattering coefficient for each 
20 of the wavelengths from the estimate; 

(iv) comparing the observed and calculated isotropic scattering 
coefficients to establish an error; 

(v) changing the estimate; and 

(vi) repeating said calculating, comparing, and changing until the 



5 



22. The method of claim 19, wherein said providing includes: 



is 



(i) 



determining an observed isotropic scattering coefficient for the 
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error reaches a desired minimum. 

23. The method of claim 22, wherein said providing includes 
selecting the estimate to generally maintain mass balance of the particles, 
5 and the estimate corresponds to an expected form of the size distribution of 
the particles. 



24. The method of claim 22, wherein the calculated isotropic 
scattering coefficient is determined with an equation having an estimated 

io product of the size distribution and. the volume fraction as an argument 

25. The method of claim 19, wherein the particle interaction 
parameter is determined from a P-Y structure factor model. 

15 26. A method of particle analysis, comprising: 

(a) exposing a number of particles to a number of light wavelengths of 
predetermined time-varying intensity; 

(b) detecting multiply scattered light from the particles in response to 
said exposing to determine a number of values each corresponding to a 

20 different one of the wavelengths, the values each being representative of a 
time-of-flight characteristic of the particles; and 

(c) providing an output determined from the values, the output 
corresponding to at least one of a particle size distribution, particle volume 
fraction, or a particle interaction parameter, the particle interaction parameter 
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corresponding to a nonlinear relationship between particle concentration and 
an isotropic scattering coefficient for the particles. 

27. The method of claim 26, wherein the particles are suspended in 
5 a fluid medium and further comprising controlling a process in accordance 

with the output. 

28. The method of claim 26, wherein the particles have a 
concentration in the fluid medium of at least about 10% by volume. 

10 

29. The method of claim 26, wherein said determining includes 
calculating at least one of volume fraction and size distribution of the particles 
in the fluid. 

15 30. The method of claim 26, wherein the particle interaction 

parameter is determined from the Percus-Yevick hard sphere model. 

31. The method of claim 30, wherein said calculating includes 
adjusting the hard sphere model to account for forces between the particles. 

20 

32. A method of analysis, comprising: ~" 

(a) exposing a fluid to an incident light, the fluid having a number of 
suspended particles therein, the suspended particles being sufficiently 
concentrated in the fluid to scatter light; 
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(b) detecting multiply scattered light in response to said exposing to 
determine a time-based value characteristic of propagation time of the 
multiply scattered light through the fluid; 

(c) determining a quantity as a function of the value, the quantity 
corresponding to an isotropic scattering coefficient; and 

(d) providing an output determined from the quantity, the output 
corresponding to at least one of particle volume fraction, particle size 
distribution, or a particle interaction parameter, the particle interaction 
parameter corresponding to particle-to-particle interactions influencing light 
scattering behavior of the particles. 

33. The method of claim 32, further comprising controlling a process 
in which the particles are altered by utilizing the output as a feedback 
variable. 

34. The method of claim 32, wherein said providing includes 
establishing an estimate corresponding to volume fraction and the size 
distribution and determining the particle interaction parameter as a function of 
the estimate. 

35. The method of claim 34, further comprising constraining the 
estimate to maintain mass balance. 

36. The method of claim 32, wherein said providing includes 
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calculating a number corresponding to the isotropic scattering coefficient as a 
function of an estimate of the volume fraction and the size distribution. 

37. The method of claim 32, wherein the incident light is intensity 
5 modulated at a predetermined frequency, and said determining includes 
comparing the incident light and the scattered light to determine the value, 
and the value is representative of a relative phase or amplitude of the 
scattered light for the predetermined frequency. 

io 38. The method of claim 37, wherein said detecting includes 

detecting the scattered light with a second sensor spaced apart from the first 
sensor by a separation distance and said determining includes calculating the 
quantity in accordance with the separation distance. 

15 39. The method of claim 37, wherein said providing includes 

determining the volume fraction or size distribution in accordance with the 
diffusion equation for multiply scattered light. 

40. The method of claim 32, wherein the particle interaction 
20 parameter is determined from a P-Y structure factor model. 

41. A system for analyzing a number of particles suspended in a 
medium in sufficient concentration to multiply scatter light, comprising: 

a light source configured to expose the medium to a number of 
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different predetermined wavelengths of incident light each having a 
predetermined time-varying intensity; 

a first sensor spaced apart from said source, said first sensor being 
configured to provide a first detection signal corresponding to detected light, 
5 the detected light being multiply scattered by the particles; 

a processor responsive to said first detection signal and being 
configured to receive an exposure signal corresponding to said incident light, 
said processor being configured to generate: (a) a number of propagation 
signals by comparing said first detection signal and said exposure signal for 
- 10 each of said wavelengths, said propagation signals each characterizing time 
of flight of the detected light through the medium resulting from multiple 
scattering by the particles for a corresponding one of said wavelengths, (b) a 
number of scattering signals each corresponding to an isotropic scattering 
coefficient of the medium and being determined from a corresponding one of 
15 said propagation signals, and (c) an output signal indicative of at least one of 
size distribution or volume fraction of the particles, said output signal being 
determined frorcusaid scattering signals and a structure factor, said structure 
factor accounting for particle-to-particle interactions influencing light scattering 
behavior of the particles; and 
20 an output device responsive to said output signal to provide an output 

corresponding to said size distribution or said volume fraction of the particles. 

42. The system of claim 41 , wherein said processor is further 
configured to determine said size distribution or said volume fraction as a 
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function of mass of the particles. 

43. The system of claim 41 1 further comprising a second sensor 
providing a second detection signal, said scattering signals being determined 
as a function of said first and second detection signals for each of the 
wavelengths. 

44. The system of claim 41 , wherein said structure factor 
corresponds to a P-Y hard sphere structure factor model. 

45. The system of claim 41 , wherein said processor further 
generates an absorption signal corresponding to an absorption coefficient of 
the medium and determines said absorption signal from said propagation 
signal. 

46. A method of particle analysis, comprising: 

(a) exposing a number of particles to a number of light wavelengths 
of predetermined time-varying intensity; 

(b) detecting multiply scattered light from the particles in response 
to said exposing to determine a number of time-based propagation 
characteristics of the particles each corresponding to a different one of the 
wavelengths; and 

(c) calculating an observed isotropic scattering coefficient for each 
of the wavelengths from the characteristics; 
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(d) determining a calculated isotropic scattering coefficient for each 
of the wavelengths from an estimate of at least one of particle size distribution 
or particle volume fraction; 

(e) comparing the observed isotropic scattering coefficient and 
5 calculated isotropic scattering coefficient for each of the wavelengths to 

establish and error; 

(f) adjusting the estimate and repeating said determining and said 
comparing until the error reaches a desired minimum; and 



10 size distribution, particle volume fraction, or a particle interaction parameter. 

47. The method of claim 46, wherein said determining includes 
establishing the particle interaction parameter as a function of the particle 



48. The method of claim 47, wherein the particle interaction 
parameter corresponds to the P-Y hard sphere structure factor. 

49. The method of claim 48, wherein the calculated isotropic 
20 scattering coefficient is determined with an equation having the structure 

factor and an estimated product of the size distribution and the volume 
fraction as arguments. 



(g) 



providing an output corresponding to at least one of the particle 



volume fraction and the particle size distribution. 



15 



50. The method of claim 46, wherein the particles include liquid 
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droplets dispersed in a fluid medium, the droplets and the medium having 
different indices of refraction. 
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where it is desired to minimize x so that it is less than or equal To an 
established error amount. Conditional 546 tests whether this minimization 
has taken place. If the minimum has been reached, then control flows to 
stage 550 and the estfmate is provided as the measured product +f(x). 
5 Subsequently, f(x) and + may be separated by performing a summation 

approximation of / $f(x)dx = $ . 
o 

However, if conditional 546 is not satisfied, then control flows to stage 
548 to improve the estimate of <J>f(x). One embodiment for rapidly converging 
to the desired minimum error is through application of a Jacobian matrix 
io relationship. Specifically, the product <K(x) is divided into N number of bins 
representing +f(x) h , where h = 1 to N. In order to update $f(x) b in each bin. 
the Jacobian matrix which describes the sensitivity of the isotropic scattering 
coefficient to changes in <t>f(x) h in each bin h and at each wavelength X, is 
computed. This Jacobian matrix is given by: 

3 = j aoKCx)), 5(<j>f(x)) l • " a(<j>f(x)) N j 
s Lac«*» rcx:» l a(*f(x» 2 ' " sotfWXj (10) 

and each element is computed through a summation approximation of 
equation (3) with *f(x) h and 4>f(x) h + A<f»f(x)„ . The updates. 6$r(x) h . 
resulting from the differences between measured and computed isotropic 
scattering coefficients, can then be obtained from: 
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.isstur iny an expected form of the PSD Generally, particulate processes can 
usually be characterized by Gaussian or log-normal distributions, both of 
which can be represented by proper parametric choices for a Weibull 
distribution described as follows: 



:(x) --- - 



cf (x-a) r • r r x - 

bL— j - > (15) 



where a. b and c describe the peak location, width, and shape, respectively, 
of the distribution. For this embodiment, d corresponds to <J> . and the 
equation (3) summation approximation becomes: 



Ac-..) 



-Ml h 



(16) 



where the expected range of particle diameter x from x min tq x max is determined 
and incremented into k number of discrete sizes x, = x min + (t-1) Ax. 
Ax = (\ ,,.,< - \ MIII1 ) / (z - t) for i=1 to z. and the terms Q Ml (\.n.x,) and g(A. r n.x, ) 
are determined from look-up tables based on k )t n, x,. 

Using Newton's method, the updates based estimation of the a t b, c,d 
parameters can be obtained from the resulting system of equations: 

3'3A4= 3 r iui;r-oo c i 

where me Jacobian matrix 3 now represents the sensitivity of isotropic 
scMtterr-g coefficients measured at wavelengths j=1 .M upon the four 
p^rsme'ers (a. b, c and d) 
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a- 5d J 



The elements of the Jacobian are numerically computed using 
equation (3). The vectors. (us')° and (Ms') c contain the 'observed and 
computed values of isotropic scattering coefficients at the wavelengths of 
interest: 

(n;>° = «>o;.g*;>:,.-<hX>, w -((h;);,(K)v-(K)L) 

(19) 

A£ is the vector updating the four parameters, a, b. c and d: 
A£ = (6a.5b.6c.5d) T (20) 
Thus, the particle sizing task now becomes to recover the four parameters (a. 
b. c and d) to describe f(x) and $ using the approach described in connection 
with loop 540. 

In still another embodiment, a low pass filter is employed to smooth the 
estimate based on the parameters for an expected form of distribution, such 
15 as the Weibull distribution. In this embodiment, the size distribution is 

computed from the parameters and then subjected to a digital low pass filter 
by averaging over a window of N"Ax. 



where £ is a factor between 0 and 1. and the summation is over the N 



'■ p "T (21) 
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intensity In other embodiments of the present invention, the light source 
arrangement is configured to provide wavelengths and time-varying intensity 
suitable for interrogation of the particular process being monitored using 
techniques known to those skilled in the art. 
5 Analyzer 150 also includes optical fiber 134 to detect re-emitted light 

from the particles in vessel 160 after multiple scattering. It is preferred that 
analyzer 150 be of rugged design suitable for industrial applications. 
Analyzer 1 50 includes a processor configured to execute software to 
determine values representative of PSD or volume fraction in accordance 
io with the present invention. Furthermore, analyzer 150 is also operatively 
coupled to control elements 170 which are used to regulate processing in 
vessel 160. By way of non-limiting example, elements 170 may include 
valves, heaters, or agitation devices electronically regulated by analyzer 150. 
Analyzer 150 is further configured with appropriate programming and 
is interfaces to provide one or more output signals to control elements 170 as a - 
function of PSD or volume fraction. Thus, system 120 provides a closed loop 
feedback control system capability for regulating chemical processes from on- 
line measurement of PSD and/or volume fraction. 



20 corresponding light source fibers 30, 130 and detecting fibers 36, 136 

geometrically arranged to approximate infinite boundary conditions such 
that light travelling from fibers 30, 130 to fibers 36, 136 does not encounter a 
substantial boundary. Equations (6) or (7) are applied to determine the 
isotropic scattering coefficient p s ' and the absorption p 4 under such infinite 



Systems 20 and 120 of FIGS. 1 and 3, respectively, are described with 
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CLAIMS TECHNOLOGY CENTER 2800 

SPECIAL PROGRAM CENTS* 

What is claimed is: 

1. A method of particle analysis, comprising: 

(a) exposing a number of particles in suspension to an incident light 
5 with a predetermined time-varying intensity, the particles being sufficiently 

close to one another to multiply scatter light; 

(b) detecting multiply scattered light from the particles in response to 
the incident light to determine a first value corresponding to an observed 
isotropic scattering coefficient of the particles: 

io (c) establishing an estimate corresponding to at least one of volume 

fraction or size distribution of the particles; 

(d) calculating a second value from the estimate, the second value 
corresponding to a calculated isotropic scattering coefficient; 

(e) comparing the first and second values to establish an error; 
15 (0 changing the estimate; and 

(9) repeating said calculating, comparing, and changing until the error 
reaches a desired minimum. 



2 The method of claim 1, wherein the incident light includes a 
jo number of different wavelengths of light, the first value is determined for each 
of the wavelengths, and said establishing, said calculating, said comparing, 
and said changing are performed for each of the different wavelengths. 
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3. The method of claim 1 , further comprising controlling a process 
as a function of the estimate and determining the estimate to generally 
maintain mass balance. 

s 

4. The method of claim 1 . wherein the estimate corresponds to 
one of a number of parameters for an expected form of the size distribution. 

5. The method of claim 4. wherein the incident light includes a 

,o number of different wavelengths of light, the first value is determined for each 
of the wavelengths, the parameters correspond to a Weibull distribution, and 

said calculating includes: 

(i) establishing a number of particle size increments; 

(ii) determining a scattering efficiency and a mean cosine of 

, s scattering angle as a function of the wavelengths and the increments; and 

(hi) performing a first summation over a range of the particle sizes 
for each of the wavelengths, the first summation having the parameters, the 
scattering efficiency, and the mean cosine of scattering angle as arguments 



20 



6 The method of claim 5. wherein said comparing includes 
perfor.n.ng a second summation over a range of the wavelengths, the second 
nation mcluding a numerical difference between the first and second 
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values as an argument. 

7. The method of claim 5, wherein said changing includes updating 
the parameters as a function of the wavelengths, the firsthand second values, 
and a third value corresponding to mass of the particles. 

8. A system for analyzing a number of particles suspended in a 
medium in sufficient concentration to multiply scatter light comprising: 

a light source with a predetermined time-varying intensity configured to 
expose the medium to a number of different predetermined wavelengths of 
incident light; 

a first sensor spaced apart from said source, said first sensor being 
configured to provide a first detection signal corresponding to detected light, 
the detected light being multiply scattered by the particles; 

a processor responsive to the first detection signal and being 
configured to receive an exposure signal corresponding to said incident light, 
said processor being configured to generate: (a) a number of comparison 
signals each corresponding to a difference with respect to time between the 
detected light and the incident light for a corresponding one of the 
wavelengths, (b) a number of scattering signals each correspondingly 
determined from the comparison signals and each corresponding to an 
oosen/ed isotropic scattering coefficient of the medium for a different one of 
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the wavelengths, and (c) an output signal indicative of one of size distribution 
or volume fraction of the particles, said output signal being determined as a 
function of said scattering signals; and 

an output device responsive to said output signal to-provide an output 
5 corresponding to the size distribution or volume fraction of the particles. 

9. The system of claim 8. further comprising a reaction vessel 
containing the particles and medium, and wherein said output device 
includes a control element operatively coupled to said reaction vessel and 

10 responsive to said output to regulate a process involving the particles. 

10. The system of claim 8. wherein said processor is further 
configured to determine said size distribution or volume fraction as a function 
of mass of the particles. 

15 

1 1. The system of claim 8, further comprising a second sensor 
providing a second detection signal, said comparison signals being 
determined as a function of said first and second detection signals for each of 
the wavelengths. 



20 



1? The system of claim 8, wherein said processor further 
determines a structure factor indicative of particle-to-particle interactions. 
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said structure factor varying in accordance with concentration of the particles in 
the medium. 

5 

13. The system of claim 8, wherein said processor is configured to 
generate a particle interaction signal representative of particle-to-particle 
interactions that vary with particle concentration and influence light scattering 
behavior of the particles. 

10 

J 

^ 14. A system for analyzing a number of particles suspended in a 

medium in sufficient concentration to multiply scatter light, comprising: 

(a) a light source with a predetermined time-varying intensity configured 
to expose the medium to a number of different wavelengths of light; 
I s (b) a sensor spaced apart from the source, said sensor being 

configured to provide a detected light signal corresponding to multiply scattered 
light from the particles at the different wavelengths in response to exposure to 
said source; 

J 

(c) a processor responsive to said detected light signal, said processor 
20 including a calculation means for generating an output signal corresponding to 
at least one of particle size distribution or volume fraction in accordance with an 
observed isotropic scattering coefficient of the medium determined from said 
detected light signal, said calculation means including a means for iteratively 
determining a structure factor from an estimate corresponding to at least one of 
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said particle size distribution or said volume fraction, said structure factor being 
representative of particle interactions that influence light scattering behavior of 
said particles; and 

5 (d) an output device responsive to said output signal to provide an 

output corresponding to at least one of said size distribution or said volume 
fraction. 

15. The system of claim 14. further comprising a reaction vessel 
10 containing the particles and the medium, and wherein said output device 

includes a control element operatively coupled to said reaction vessel and 
responsive to said output to regulate a reaction involving the particles. 

16. The system of claim 14 f wherein said calculation means includes 
1 5 an estimating means for iteratively determining said size distribution or said 

volume fraction as a function of mass of the particles. 



1 7. The system of claim 14, wherein said structure factor is 
dependent on said particle size distribution and said particle volume fraction. 

20 and corresponds to a P-Y hard sphere model. 

18. The system of claim 14. wherein said calculation means includes 
a number of parameters corresponding to a Weibull distribution. 



AMENDED SHEET 



r 



SENT BY: WOODARD EMHARDT ; 




1-22-99 2:56PM; 317 634 3456 ^ 




7033057724; #14/2*: 



PCT/US9V /20 53S 
IPEA/US 2 2 JAN 1999 



64 



19. A method of particle analysis, comprising: 

(a) exposing a medium with a number of suspended particles to a 
number of light wavelengths, the wavelengths each being intensity-modulated 

5 at a predetermined frequency; 

(b) detecting multiply scattered light from the medium in response to 
said exposing to characterize propagation of the multiply scattered light 
through the medium with a number of values, the values each corresponding 
to a different one of the wavelengths and each being representative of at 

io least one of a phase or an amplitude of the multiply scattered light relative to 
the predetermined frequency; and 

(c) providing an output determined from the values, the output 
corresponding to at least one of a particle size distribution, particle volume 
fraction, or a particle interaction parameter, the particle interaction parameter 

is corresponding to a nonlinear relationship between particle concentration and 
an isotropic scattering coefficient for the particles. 



in accordance with the output. 

20 

21. The method of claim 19, wherein said providing includes 
determining the observed isotropic scattering coefficient and an absorption 



20. 



The method of claim 19, further comprising controlling a process 
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coefficient for the particles for each of the wavelengths from a corresponding 
one of the values. 

22. The method of claim 19. wherein said providing includes: 
5 (i) determining an observed isotropic scattering coefficient for the 

particles for each of the wavelengths from a corresponding one of the values; 

(ii) establishing an estimate corresponding to at least one of the 

j size distribution or the volume fraction; 

(iii) determining a calculated isotropic scattering coefficient for each 

io of the wavelengths from the estimate; 

(iv) comparing the observed and calculated isotropic scattering 

coefficients to establish an error: 

(v) changing the estimate; and 

(vi) repeating said calculating, comparing, and changing until the 
15 error reaches a desired minimum. 



23. The method of claim 22, wherein said providing includes 
selecting the estimate to generally maintain mass balance of the particles, 
and the estimate corresponds to an expected form of the size distribution of 
the particles 

»A. The method of claim 22. wherein the calculated isotropic 
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scattering coefficient is determined with an equation having an estimated 
product of the size distribution and the volume fraction as an argument. 

5 25. The method of claim 19, wherein the particle interaction 

parameter is determined from a P-Y structure factor model: 

26. A method of particle analysis, comprising: 

(a) exposing a number of particles to a number of light wavelengths of 
1 0 predetermined time-varying intensity; 

(b) detecting multiply scattered light from the particles in response to 
said exposing to determine a number of values each corresponding to a 
different one of the wavelengths, the values each being representative of a 
time-of-flight characteristic of the particles; and 

15 (c) providing an output determined from the values, the output 

corresponding to at least one of a particle size distribution or volume fraction 
and being determined in accordance with a particle interaction parameter, the 
particle interaction parameter being representative of a nonlinear relationship 
between particle concentration and an isotropic scattering coefficient for the 



27. The method of claim 26, wherein the particles are suspended in 
a fluid medium and further comprising controlling a process in accordance with 
the output. 



20 



particles. 
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28. The method of claim 26. wherein the particles have a 
concentration in the fluid medium of at least about 10% by volume. 

29. The method of claim 26. wherein said determining includes 
calculating at least one of volume fraction and size distribution of the particles 
in the fluid. 

30. The method of claim 26. wherein the particle interaction 
parameter is determined from the Percus-Yevick hard sphere model. 

31 . The method of claim 30. wherein said calculating includes 
adjusting the hard sphere model to account for forces between the particles. 



l5 32. A method of analysis, comprising: 

(a) exposing a fluid to an incident light, the fluid having a number of 
suspended particles therein, the suspended particles being sufficiently 
concentrated in the fluid to scatter light; 

(b) detecting multiply scattered light in response to said exposing to 
20 determine a time-based value characteristic of propagation time of the multiply 

scattered light through the fluid. 

(c) determining a quantity as a function of the value, the quantity 
corresponding to an isotropic scattering coefficient; and 

(d) providing an output corresponding to at least one of volume 
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fraction, particle size distribution, or a particle interaction parameter, the particle 
interaction parameter corresponding to particle-to*particle interactions 
influencing light scattering behavior of the particles, said providing including 
5 calculating a number representative of the isotropic scattering coefficient as a 
function of an estimate corresponding to the volume fraction and the size 
distribution. 

33. The method of claim 32. further comprising controlling a process 

^ 10 in which the particles are altered by utilizing the output as a feedback variable. 

* ■? 

34. The method of claim 32, wherein said providing includes 
establishing an estimate corresponding to volume fraction and the size 
distribution and determining the particle interaction parameter as a function of 

15 the estimate. 

35. The method of claim 34 p further comprising constraining the 
, estimate to maintain mass balance. 
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37. The method of claim 32, wherein the incident light is intensity 
modulated at a predetermined frequency, and said determining includes 
comparing the incident light and the scattered light to determine the value. 

s and the value is representative of a relative phase or amplitude of the 
scattered light for the predetermined frequency. 

38. The method of claim 37, wherein said detecting includes 
detecting the scattered light with a second sensor spaced apart from the first 

io sensor by a separation distance and said determining includes calculating the 
quantity in accordance with the separation distance. 

39. The method of claim 37 f wherein said providing includes 
determining the volume fraction or size distribution in accordance with the 

15 diffusion equation for multiply scattered light. 

40. The method of claim 32, wherein the particle interaction 
parameter is determined from a P-Y structure factor model. 

20 4 1. A system for analyzing a number of particles suspended in a 

medium in sufficient concentration to multiply scatter light, comprising: 
a light source configured to expose the medium to a number of 
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different predetermined wavelengths of incident light each having a 
predetermined time-varying intensity; 

a first sensor spaced apart from said source, said first sensor being 
configured to provide a first detection signal corresponding*4o detected light. 
5 the detected fight being multiply scattered by the particles; 

a processor responsive to said first detection signal and being 
configured to receive an exposure signal corresponding to said incident light. 
^ said processor being configured to generate: (a) a number of propagation 

signals by comparing said first detection signal and said exposure signal for 
io each of said wavelengths, said propagation signals each characterizing time 
of flight of the detected light through the medium resulting from multiple 
scattering by the particles for a corresponding one of said wavelengths, (b) a 
number of scattering signals each corresponding to an isotropic scattering 
coefficient of the medium and being determined from a corresponding one of 
is said propagation signals, and (c) an output signal indicative of at least one of 

i 

^ size distribution or volume fraction of the particles, said output signal being 

determined from said scattering signals and a structure factor, said structure 
factor accounting for particle-to-particle interactions influencing light scattering 
behavior of the particles; and 
10 an output device responsive to said output signal to provide an output 

corresponding to said size distribution or said volume fraction of the particles. 
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42. The system of claim 41, wherein said processor is further 
configured to determine said size distribution or said volume fraction as a 
function of mass of the particles. 

43. The system of claim 41, further comprising a second sensor 
providing a second detection signal, said scattering signals being determined 
as a function of said first and second detection signals for each of the 
wavelengths. 



44. The system of claim 41 , wherein said structure factor 
corresponds to a P-Y hard sphere structure factor model. 

45. The system of claim 41. wherein said processor further 
generates an absorption signal corresponding to an absorption coefficient of 
the medium and determines said absorption signal from said propagation 

signal. 



46. A method of particle analysis, comprising: 

(a) exposing a number of particles to a number of light wavelengths 
of predetermined time-varying intensity; 

(b) detecting multiply scattered light from the particles in response 
to said exposing to determine a number of time-based propagation 
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characteristics of the particles each corresponding to a different one of the 

wavelengths; and 

( c) calculating an observed isotropic scattering coefficient for each 

of the wavelengths from the characteristics; 
5 (d ) determining a calculated isotropic scattering coefficient for each 

of the wavelengths from an estimate of at least one of particle size distribution 

or particle volume fraction; 

(e) comparing the observed isotropic scattering coefficient and 
calculated isotropic scattering coefficient for each of the wavelengths to 

io establish and error; 

(f, adjusting the estimate and repeating said determining and said 

comparing until the error reaches a desired minimum; and 

(g) providing an output corresponding to at least one of the particle 
size d.stribution. particle volume fraction, or a particle interaction parameter. 



I 5 

J 
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47. The method of claim 46. wherein said determining includes 
establishing the particle interaction parameter as a function of the particle 
volume fraction and the particle size distribution. 

48. The method of claim 47. wherein the particle interaction 
parameter corresponds to the P-Y hard sphere structure f-ctor. 
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49. The method of cfaim 48. wherein the calculated isotropic 
scattering coefficient is determined with an equation having the structure 
factor and an estimated product of the size distribution and the volume 
fraction as arguments. 

50. The method of claim 46. wherein the particles include liquid 
droplets dispersed in a fluid medium, the droplets and the medium having 

""^ different indices of refraction. 
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